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Abstract 

Local SU(5) F-theory models lead naturally to Yukawa couplings for the third gen- 
eration of quarks and leptons, but inducing Yukawas for the lighter generations has 
proven elusive. Corrections coming from gauge fluxes fail to generate the required 
couplings, and naively the same applies to instanton effects. We nevertheless revisit 
the effect of instantons in F-theory GUT constructions and find that contributions 
previously ignored in the literature induce the leading non-perturbative corrections 
to the Yukawa couplings. We apply our results to the case of 10 x 5 x 5 couplings 
in local £77(5) F-theory GUTs, showing that non-perturbative effects naturally lead 
to hierarchical Yukawas. The hypercharge flux required to break SU(5) down to 
the SM does not affect the holomorphic Yukawas but does modify the profile of the 
wavef unctions, explaining the difference between the D-quark and lepton couplings 
at the unification scale. The combination of non-perturbative corrections and mag- 
netic fluxes allows to describe the measured lepton and D-quark masses of second 
and third generations in the SM. 
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1 Introduction 

One of the most difficult puzzles of the Standard Model (SM) is the structure of fermion 
masses and mixings. If the SM arises as a low-energy limit of an underlying string the- 
ory [1] , it should be possible to understand this structure in terms of more fundamental 
parameters characterizing the string vacuum. In particular the values of Yukawa cou- 
plings in string compactifications are determined by the geometric properties of extra 
compactified dimensions. An explicit computation of Yukawa coupling constants seems 
then quite difficult since we would need detailed information about the geometry of the 
corresponding compact space. 

It has been however realized [2] that, due to the localization properties of branes, some 
quantities of physical interest do not depend on the full geometry of the compactification 
space but rather on local information around the region in which the SM fields are local- 
ized. This is particularly the case of Type IIB string compactifications with the SM fields 
localized on Dp-branes with p < 7 and also on F-theory constructions. In these Type 
IIB compactifications the Yukawa couplings are obtained as overlap integrals involving 
the three wavefunctions of the quark/lepton and Higgs states involved. In particular, 
in the case of F-theory SU(5) GUT's the quark/lepton fields are localized on complex 
matter curves with 5, 10 quantum numbers. Yukawa couplings appear at the points of 
intersection of three matter curves corresponding to a right-handed fermion, a left-handed 
fermion and a Higgs multiplet transforming also as a 5-plet. With the wavefunctions of 
the three fields localized on the matter curves, the overlap integral is then dominated by 
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the local properties of the wavefunction around the intersection point and hence no global 
information is required to compute the holomorphic Yukawa couplings. This opens the 
door to the explicit computation of Yukawa couplings in string compactifications with 
non-trivial curved compact spaces. 

Particularly interesting from this point of view are the mentioned local F-theory £77(5) 
GUT models [3-6], which have been recently the subject of intense study (for reviews see 
e.g. [7]). These models are able to combine advantages of heterotic compactifications 
(gauge coupling unification) and Type IIB orientifolds (localization of the SM fields and 
moduli fixing through closed string fluxes). In these F-theory models the local dynamics 
on the matter curves is governed by the 8d effective action of 7-branes, and one can obtain 
explicit local expressions for the wavefunctions of the matter fields. The Yukawa couplings 
arise at the triple intersection of matter curves [3-6,8,9] , and can be computed from a 
superpotential of the form W ~ J Tr (F A $), where F = dA—iAAA is the field strength of 
the 8-dimensional gauge fields and $ is a field parametrizing fluctuations in the transverse 
dimensions to the branes, and the integral extends over the 7-brane worldvolume where 
the £77(5) degrees of freedom are localized. Within this simple scheme one finds that only 
one generation of quark/lepton fields gets a Yukawa coupling and may eventually become 
massive [10,11]. This is analogous to the result obtained in Type II toroidal orientifolds 
in which the Yukawas may be computed explicitly [15-17]. This is an interesting starting 
point since indeed in the SM the third generation is much heavier than the rest and one 
may think that some additional corrections could give masses to the first two generations. 1 
It was first thought [10] that the presence of the world- volume fluxes required both to get 
chirality and break the SU(5) symmetry down to the SM could be the source of these 
corrections. However it was soon realized [18-20] that open string fluxes do not modify 
at all the holomorphic Yukawa couplings and hence cannot give rise to Yukawa couplings 
for the lighter generations. 

In [21,22] it was pointed out that non-perturbative effects from distant D3-instantons 

could be the source of the required corrections. The most obvious such corrections where 

found to be proportional to e J Tr ($F A F) and turn out to have an alternative useful 

1 For different approaches to the generation of hierarchies of fermion masses in F-theory unification 
see e.g. [12-14]. 
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description in terms of non-commutative geometry [18]. Indeed such corrections where 
shown to lead to the required corrections in a simple toy model with gauge group U(3) 
[22]. It was however already pointed out in this reference that in realistic cases, namely 
£77(5) GUT's with an enhanced symmetry group £0(12) or Eq j7j8 at the Yukawa triple 
intersection point, such corrections identically vanish, since the cubic trace is zero in 
orthogonal and exceptional groups. This seemed again to make problematic the generation 
of fermion hierarchies in F-theory GUT constructions. 

In this paper we reexamine all these issues and point out that non-perturbative D3- 
instanton effects give rise to additional corrections, some of them previously overlooked. 
Such corrections to the superpotential have the form 



where e is a small parameter and 9 n are holomorphic functions of the local coordinates. 
As mentioned above the contribution with n — 1 vanishes for the realistic cases where the 
Yukawa enhancement groups are £0(12) or Eq. We hence study in detail the remaining 
leading corrections to the Yukawa couplings induced by the n = and n = 2 terms, 
applied to the £0(12) case which is relevant for the Yukawa couplings of charged leptons 
and D-quarks. 

Describing non-perturbative corrections as in (1.1) simplifies the procedure to compute 
corrected Yukawas. More precisely, one may apply dimensional reduction techniques to 
express them in terms of a triple overlap of zero mode internal wavefunctions. In this 
sense, the presence of W np has a two-fold effect. On the one hand it modifies the zero 
mode internal wavefunction profile and on the other hand it induces new 8d couplings that 
upon dimensional reduction become new sources of Yukawa couplings. As in [22], taking 
both effects into account gives an interesting Yukawa pattern, in which the holomorphic 
Yukawas depend on 9 n but are independent of worldvolume fluxes. While in [22] this result 
can be guessed based on a dual non-commutative description of the 7-brane superpotential 
[18,21], for the general case (1.1) such description is not available. Nevertheless, one can 
still generalize the results of [18] to obtain a residue formula that computes holomorphic 
Yukawas, and where the flux-independence of the latter is manifest. 

Interestingly enough we find that a hierarchy of mass eigenvalues of the form (1, e, e 2 ) 




(1.1) 
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is automatically present, explaining the observed hierarchical structure. Here e ~ 10~ 2 
is a small non-perturbative parameter measuring the size of the effects induced by the 
distant instantons. The holomorphic Yukawa couplings obtained are identical for both D- 
quarks and leptons, since they live in the same SU (5) representations and the holomorphic 
Yukawas are flux independent. This looks problematic since running up in energies the 
observed D-quark and lepton masses, unification of Yukawa couplings Y b Sjd = Y T ^ e does 
not hold experimentally, rather leptons of the second and third generations tend to have 
larger Yukawa couplings than the respective D-quarks at the unification scale. We find 
however that the hypercharge flux required for the SU (5) — > SM symmetry breaking may 
explain this difference. Roughly speaking, the difference may be understood as arising 
from the fact that the wavefunctions for leptons are more localized than those of D-quarks, 
due to the fact that they have larger hypercharge quantum numbers. 

The structure of this paper is as follows. In Section 2 we review the construction 
of local F-theory GUTs. In section 3 we construct a local SU (5) GUT model with en- 
hanced SO(12) symmetry, which describes the Yukawa couplings of charged leptons and 
D-quarks. The spectrum of zero modes reproduces the matter content of the MSSM, but 
the Yukawa couplings exhibit the rank-one structure mentioned above. In section 4 we 
introduce the non-perturbative effects that will give rise to the superpotential (1.1), and 
compute the corrected zero mode equations. Such equations are solved in section 5 for 
the model of SO (12) enhancement constructed before, while the corresponding Yukawas 
are computed in section 6. The discussion of these last two sections is slightly techni- 
cal, and the reader not interested in such details may safely skip to section 7, where a 
phenomenological analysis of the final Yukawa couplings is performed. In particular, we 
confront our results with the measured masses of D-quarks and charged leptons, showing 
that a natural hierarchy of masses arises and that the effect of the hypercharge flux allows 
us to understand the ratios between them. 

Several technical details have been relegated to the appendices. Appendix A solves 
the zero mode wavefunctions for the SO{12) model in absence of non-perturbative effects, 
and compute the wavefunction normalization factors which encode the hypercharge flux 
dependence of the Yukawa couplings. Appendix B discusses in some detail the choice 
of worldvolume fluxes made for this 5*0(12) model, motivating them via the notion of 
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local chirality in F-theory. Appendix C derives the non-perturbative superpotential (1.1), 
and shows that the D-term is not corrected. Finally, in appendix D we derive a residue 
formula for the non-perturbative Yukawa couplings, that allows to cross-check and extend 
the results obtained in the main text. 

2 Review of local F-theory models 

Following the general scheme of [3-6] (see also [23-30]), in order to construct a local F- 
theory GUT model one may consider a stack of 7-branes wrapping a compact divisor Scut 
of the threefold base B of an elliptically-fibered Calabi-Yau fourfold. The gauge degrees 
of freedom that arise from Sqvt are specified by the particular set of (p, q) 7-branes that 
are wrapped on Scut or, in geometrical terms, by the singularity type of the elliptic fiber 
on top of such 4-cycle. Hence, one may easily engineer local models where the GUT gauge 
group Ggut is given by SU(5), 5*0(10) or even E$. 

Besides the stack of 7-branes on Scut, a semi- realistic F-theory model will contain 
further 7-branes that wrap another set of divisors Si, which intersect Scut on certain 
curves Ej. On top of the latter set of curves of Scut the singularity type of the elliptic 
fiber is enhanced, in the sense that the Dynkin diagram that is associated to the singularity 
corresponds to a higher rank Lie group that contains Ggut- In practice, this implies 
that new degrees of freedom appear at the intersection of the 7-branes, more precisely 
chiral matter multiplets in a certain representation of Ggut, localized at the so-called 
matter curves Ej. 

Finally, two or more matter curves may meet at a point p G Scut and at that point 
the singularity is promoted to an even higher one, such that the corresponding Lie group 
G p not only contains Gqut but also each of the G^ involved. This time there are no new 
degrees of freedom arising at the point p, but rather contact interactions involving the 
chiral multiplets from each curve Ej. Of particular interest are those cases where three 
matter curves meet at p, as they give rise to Yukawa couplings between chiral multiplets 
of the GUT matter fields. 

Of course, in the process of describing a local model one must not only specify the 
gauge group Gqvt, but also the enhanced group G^ at each of the matter curves. This 
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information and the intersection loci of matter curves determines the groups G p at each 
point p where Yukawa couplings develop. Typically, starting with a GUT gauge group 
such as Ggut = SU (5) one may end up with enhanced groups at Yukawa points p such as 
SO (12), Eq, E 7 or E s . In the next section we will analyze a local model that describes the 
case where Ggut = SU(5) and G p = SO (12), which corresponds to the setup describing 
down-type Yukawas for a local 577(5) F-theory model. 

While the above geometric picture is already quite illuminating, one of the most power- 
ful results of [3-6] is to provide a simple framework to compute the matter content arising 
at each curve £, and the Yukawa couplings at their triple intersections. Such framework 
makes use of a 8d effective action related to a stack of 7-branes which, upon dimensional 
reduction on a 4-cycle S, provides all the dynamics of the 4d degrees of freedom [4]. 2 In 
particular, the Yukawa couplings between 4d chiral fields arise from the superpotential 



where m* is the F-theory characteristic scale, F = dA — %A A A is the field strength of 
the 8d gauge vector boson A arising from 7-branes, and $ is a (2,0)-form on the 4-cycle 
S describing its transverse geometrical deformations. Near the Yukawa point p, we can 
take A and $ to transform in the adjoint of the enhanced non-Abelian group G p , which 
in our case will be given by SO(12). Further dynamics of this system is encoded in the 
D-term 



where u stands for the fundamental form of S. Together with the superpotential, this D- 
term relates the spectrum of 4d zero and massive modes to a set of internal wavefunctions 
along S, and the couplings between these 4d modes to the overlapping integrals of such 
wavefunctions. 

Notice that from this latter viewpoint we seem to have a single divisor S with a higher 
gauge group G p . One must however take into account that both $ and A have a non- 
trivial profile. On the one hand the nontrivial profile for $ (more precisely the fact that 
the rank of ($) jumps at the curves £j) takes into account the fact that we do not have 
a single divisor S, but rather a set of intersecting divisors Scut and Si. A non- vanishing 
2 See [19] for a derivation in terms of a 8d SYM Lagrangian. 




(2.1) 




(2.2) 
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(<&) then breaks the would be gauge group G p to the subgroup Ggut x IL w hh Gj 
the gauge groups of the 7-branes wrapping the divisors 5j, typically chosen to be U(l). 

On the other hand, the effect of (A) is to provide a 4d chiral spectrum and to further 
break the GUT gauge group G G ut down to the subgroup that commutes with (A), as 
it is usual in compactifications with magnetized D-branes [16,31-34]. As a result, one 
may obtain a 4d MSSM spectrum from the above construction by first engineering the 
appropriate GUT 4d chiral spectrum via (<£>) and an (A) which commutes with Ggut, 
and by then turning on an extra component of (A) along the hypercharge generator Q Y 
in order to break Ggut —> Gmssm [5]. Generically, the presence of a non- vanishing field 
strength F Y along the hypercharge generator is the only way to break the GUT gauge 
group down to the MSSM one. As a result, all the physics of the MSSM that differ from 
the parent GUT physics must depend on the data that describe (Fy)- 

Finally, in addition to the above set of divisors hosting the MSSM gauge and matter 
content, there will be in general other divisors also wrapped by branes which may source 
non-perturbative effects. Typical examples are 7-branes with a gauge hidden sector of the 
theory that undergoes a gaugino condensate, or Euclidean 3-branes with the appropriate 
structure of zero modes to contribute the the superpotential of the 4d effective theory. 
Such ingredients are usually not considered in the construction of F-theory local models, 
and indeed they will not be present in the 5*0(12) model described in the next section. 
However, as we will review in section 4, they are crucial in endowing F-theory local models 
with more realistic Yukawa couplings. In fact, one of the main results of this work is to 
show this point for the class of SO(12) local models that we now proceed to describe. 

3 The SO(12) model 

In this section we describe in detail the 50(12) local model which we will analyze in the 
rest of the paper. Following the common practice in the F-theory literature, we will first 
specify the structure of 7-brane intersections and matter curves that breaks the 50(12) 
symmetry down to SU(5) x U(l) 2 , and then add the worldvolume flux that induces 4d 
chirality and breaks the SU (5) GUT spectrum down to the MSSM. 

While in this section we will use the language of F-theory local models, it is important 
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to notice that the model at hand admits a more intuitive description in the framework of 
intersecting D7-branes in type IIB orientifolds. We will exploit such vantage point in the 
next section, in order to gain some insight on the non-perturbative corrections that can 
affect our local model. 

3.1 Matter curves 

Following the general framework described in the previous section, let us consider a local 
model where the symmetry group at the intersection point of three matter curves is 
G p = SO (12). Away from this point, this group is broken to a subgroup because (<£>) 7^ 0. 
One can then engineer a (<&) such that generically SO (12) is broken to SU(5) x U(l) 2 , 
except for some complex curves where there is an enhancement to either SO (10) x U(l) 
or SU(6) x U(l). In this way, we can identify G s = SU(5) as the GUT gauge group and 
the enhancement curves as matter curves where chiral matter wavefunctions are localized. 

In order to make the above picture more precise let us consider the generators of 
5*0(12), in terms of which we can express the particle spectrum of our local GUT model. 
These generators can be decomposed as {Hi, E p }, where the H, h % — 1, • • • , 6, belong to 
the Cartan subalgebra of SO (12) and the E p are step generators. 3 Recall that 

[H i ,E„]=p i E p (3.1) 

where pi is the i-th component of the root p. The 60 non-trivial roots are given by 

( ±1, ±1,0, 0,0,0 ) (3.2) 

where the underlying means all possible permutations of the vector entries. 

Let us now choose the vev of the transverse position field $ = Q xy dx A dy to be 

($ xy ) =m 2 (xQ x + yQ y ) (3.3) 

where m 2 is related to the intersection slope between 7-branes as explained in section 

4, and the charge operators Q x and Q y are the following combinations of generators of 

3 Throughout this work we use the standard form of the SO(2N) generators in the fundamental 
representation [35]. 
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elements of the 50(12) Cartan subalgebra 



Qx 




(3.4) 



This choice of ($) describes a 50(12) local model that is similar to the U(3) toy model 
analyzed in [22] in several aspects. This will allow us to apply several useful results of [22] 
to the more realistic case at hand. 

Given (3.3) one can understand the 50(12) symmetry breaking pattern described 
above as follows. In general the step generators E p satisfy 



with g$ a holomorphic function of the complex coordinates x, y of the 4-cycle 5. The 
subgroup of 50(12) not broken by the presence of this vev corresponds to those generators 
that commute with (<E>) at any point in 5. This set is given by the Cartan subalgebra of 
50(12) and to those step generators E p such that q<s>(p) = for all x,y. It is easy to see 
that such unbroken roots are given by 



together with the Cartan generators. Therefore, from the symmetry group 50(12) only 
the subgroup SU(5) x U(l) 2 remains as a gauge symmetry, and we can identify Gs = 
SU(5) as our GUT gauge group. 

On the other hand, the broken generators of 50(12), that have q$ ^ for generic x, y, 
allow us to understand the pattern of matter curves and to classify the charged matter 
localized therein. Such broken roots and their charges g$ are displayed in table 1. 

From this table we see that there are three complex curves within 5 where the bulk 
symmetry SU(5) x U(l) 2 is enhanced, in the sense that there g$ = for an additional set 
of roots. Concretely, for x = there are 10 additional roots that together with those in 
(3.6) complete the SU(6) root system. We have labeled such matter curve as E a , so that 
in the language of the previous section we would have that G Sa = SU(6) x U(l). These 
extra set of roots whose q$ vanishes at E a can be split into subsets that have different 
g$ away from E a . It is easy to convince oneself that each of these subsectors must fall 
into complete weight representations of SU(5), which in turn correspond to the matter 



[{Qxy),E p ] = m 2 q${p)E p 



(3.5) 



(0,1,-1,0,0,0) 



(3.6) 
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Table 1 : Data of broken generators 

localized at the curve. In the case of S a , there are two sectors a + and a~ that correspond 
to the representations 5 and 5 of SU(5), respectively, as shown in table 1. 

Similarly to S a , at the curve = {y = 0} there are 20 extra unbroken roots and 
SU(5) x U{1) 2 is enhanced to 5*0(10) x U(l), giving rise to the representations 10 and 
10. The third matter curve is given by S c = {x = y}, where there is also an enhancement 
to 51/(6) x 1/(1). 

Finally, let us consider a set of quantities that only depend on each root sector of the 
model. These are the symmetrized products 4 

S mn (p) = S(E p Q m Q n ) ; m,n = x,y (3.7) 

where the generators are taken in the fundamental representation of 5*0(12). As we will 
see, the equations of motion satisfied by the zero modes at the matter curves will depend 
on these quantities. For the broken roots we obtain 

Smn(p) = S mn (p)E p (3.8) 

where the s mn = s nm are constants also displayed in table 1. 
4 S(Ai ■ ■ ■ An) = (^i ' ' ' An + all permutations) 
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3.2 Worldvolume flux 

To obtain a 4d chiral model the above pattern of matter curves is not enough, and it 
is necessary to add a non-trivial background worldvolume flux (F) to our local F-theory 
model. Just like the position field, such flux is usually chosen along the Cartan subalgebra 
of SO (12), so that it commutes with (& X y) an d the equations of motion of our system are 
simplified. Moreoever, considering a component of (F) along the hypercharge generator 
Q Y allows to break the GUT gauge group SU(5) down to SU(3) x SU(2) x U(1) Y , this 
being in fact the only way to achieve GUT symmetry breaking for the most generic class 
of F-theory GUT models. 

In order to construct a worldvolume flux with the desired properties we proceed in 
three steps. First we add a flux (Fi) analog of the one introduced in the U(3) toy model 
of [22]. Just like in there, this flux will create chirality on the curves S a and E&, selecting 
the sectors a + and b + as the ones that contain the chiral matter of the model, as opposed 
to a~ and b~ . Then we add an extra piece (F 2 ) such that the matter curve S c also 
contains a chiral spectrum: a typical requirement to achieve an acceptable Higgs sector. 
None of these previous fluxes further break the gauge group SU(5) x U(l) 2 so, finally, 
we include a flux (F Y ) along the hypercharge generator Q Y that breaks SU (5) down to 
SU(3) x SU(2) x U(1) Y . 

To proceed we then consider the flux 

(Fi) = i (M x dx A dx + M y dy A dy) Q F (3.9) 

where 

Qf = \(H± -H 2 -H 3 -H 4 -H 5 - H 6 ) = -Q x - Q y (3.10) 

which is the analog of the flux introduced in the U(3) toy model of [22]. To analyze the 
effect of this flux it is convenient to define the Q^-charge of the roots E p according to 

[Q F ,E p } = q F (p)E p (3.11) 

The roots in (3.6) are clearly neutral under this flux component (Fi), and so the gauge 
symmetry SU (5) is not broken further by its presence. The roots in the sectors a and b 
are however not neutral. Hence, if the integral of (3.9) over each of these curves does not 
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vanish, they will each host a chiral sector of the theory. In the following we will assume 
that this is the case and that induces a net chiral spectrum of three 5's in the curve 
S a and three 10's in the curve S 6 . If this chiral spectrum can be understood in terms of 
local zero modes in the sense of [36], then such chiral modes should arise in the sectors 
a + and b + of table 1, respectively, and by the results of appendix A one should choose 
M x < < My to describe them locally. 

Notice that the roots belonging to the c sector are neutral under (3.10), and so the 
spectrum arising from the curve S c is unaffected by the presence of As the SO (12) 
triple intersection point is where down-like Yukawa couplings arise from, we do need one 
5 in such curve, but however no 5 so that no undesired 55 mass terms appear. This 
chiral spectrum on the sector c can be achieved by adding the following extra piece of 
worldvolume flux 



It is easy to check that the particles localized at the matter curve S c are now non-trivially 
charged under the flux background, and that a local chiral spectrum can be achieved if we 
choose N a 7^ Nf,. In particular, as shown in appendix B for N a > iV b one obtains net local 
chirality in the sector c + , yielding the desired 5 which is the SU(5) down Higgs. Notice 
that those particles at the curves S a and S b are also charged under (3.12). However, by 
construction the number of (local) families in such curves is independent of the flux (F 2 ), 
as also shown in appendix B. 

Let us finally add a third piece of worldvolume flux which, unlike (3.9) and (3.12), will 
break the SU(5) gauge group down to the MSSM. As usual, such flux should be turned 
along the hypercharge generator Qy, and a rather general choice is given by 



We have chosen the hypercharge flux to be a primitive (l,l)-form, so that it satisfies 
automatically the equations of motion for the background. Note that (3.13) has two 
components that are easily comparable with the previous flux components (3.9) and (3.12). 



(F 2 ) = i (dx Ady + dyA dx) (N a Q x + N b Q y ) 



(3.12) 



(F Y ) = % (dx A dy + dy A dx) N Y + (dy A dy - dx A dx) N Y Qy (3.13) 



where 




(3.14) 
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The first component, proportional to the flux density Ny, is quite similar to (F 2 ). Indeed, 
as happens for (3.12) its pullback vanishes over the matter curves £ a and E&, and so it 
does not contribute to the (local) index that computes the number of chiral families in the 
sectors a and b. The second component, proportional to Ny, ma Y m principle affect the 
chiral index over the curves £ a and E& but, following the common practice in the GUT 
F-theory literature, we will assume that this is not the case. Globally one requires that 

/ (Fy) = I (F Y ) = (3.15) 

so that three complete families of quarks and leptons remain at the curves £ a and £ b 
after introducing the hypercharge flux. Locally, we demand that the local zero modes still 
arise from the sectors a + and b + , and this amounts to require flux densities such that 

M x + q Y N Y < < M y + q Y N Y (3.16) 

for every possible hypercharge value qy in the sectors a + and b + , see table 2 below. 

While innocuous for the matter spectrum at the curves E Q , £&, the hypercharge flux 
is supposed to modify the chiral spectrum of curve E c , in order to avoid the doublet- 
triplet splitting problem of SU(5) GUT models [5]. Indeed, one typically assumes that 
J c {Fy) 7^ 0, and since (3.13) couples differently to particles with different hypercharge, 
this implies a different chiral index for the doublet and for the triplet of 5. Locally, we 
have that the total flux seen near the Yukawa point by the doublets on the sector c + is 

F tot , 2 = N Y + 2(N a -N b ) (3.17) 

while the flux seen by the triplets is 

Ftot,3 = -\n y + 2(N a - N b ) (3.18) 

Hence, in order to have a vector-like sector of triplets in the local model we can set 

N Y = 3(N a - N b ) (3.19) 

and then assume that such vector-like spectrum is massive. Notice that this condition still 
yields a chiral sector for the doublets and so forbids a /x-term for them. Indeed, imposing 
(3.19) we have that 

F tot , 2 = ^N Y (3.20) 
14 



which in general will induce a net chiral spectrum of doublets in the curve E c . Hence, 
imposing (3.19) the combined effect of (F 2 ) and (Fy) is such that doublets of 5 in the 
sector c feel a net flux, while triplets do not. One may then choose the flux density N Y 
such that it yields a single pair of MSSM down Higgses at the curve E c . 

To summarize, the total worldvolume flux on this local SO (12) model is given by 

(F) = i(dy Ady - dx Adx)Q P + i(dx Ady + dy Adx)Q s (3.21) 
+ i(dy Ady + dx A dx)M xy Q F 

where we have defined 

Q P = MQ F + N Y Q Y (3.22) 
Qs = N a Q x + N b Q y + N y Qy (3.23) 

and 

M = X -(M y - M x ) M xy ee i(M„ + M x ) (3.24) 

Note that the combination of flux densities M xy corresponds to an Fl-term, which will be 
set to vanish whenever supersymmetry is imposed. 

Just like in [22], we can now express the vev of the corresponding vector potential A 
in the holomorphic gauge defined in [11], namely as 

(A) ho1 = i [(Q P - M xy Q F )x - yQs] dx - % [(Q P + M xy Q F )y + xQ s ] dy (3.25) 

this being the quantity that will enter into the equation of motion for the zero mode 
wavefunctions at the curves E a , E 6 and E c . The combined effect of the background (<&) 
and (A) breaks 5*0(12) to SU(3) x SU(2) x U(l) 3 , and as a result the sectors a, b and c 
split into further subsectors compared to table 1. The content of charged particles under 
the surviving gauge group is shown in table 2, where we have also displayed the charges 
of each sector under the operators Qp, Q x , Q y and Qy. We have also included the values 
of qs and qp, which are defined as 

[Qs, E p \ = q s (p)E p [Q P , E p \ = q P (p)E p (3.26) 

and which, unlike the other charges, depend on the flux densities of the model. As 
discussed below and in appendix A, each of these sectors obeys a different zero mode 
equation, and so it is described by a different wavefunction. 
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Sector 


Root 


517(3) 


SU(2) 


QY 


IF 


q x 


q y 


qs 


qp 


at 


(1,-1,0,0,0,0) 


3 


1 


l 

3 


1 


-1 





-N a - iAV 


M-\N Y 


a i 


(-1,1,0,0,0,0) 


3 


1 


1 
3 


-1 


1 





N a + ±N Y 


— M + |iVy 




(1,0,0,0,-1,0) 


1 


2 


1 
2 


1 


-1 





-N a + \N Y 


M+ \N Y 




(-1,0,0,0,^0) 


1 


2 


1 
2 


-1 


1 





N a - \N Y 


—M - |iVy 


K 


(0,1,1,0,0,0) 


3 


1 


2 
3 


-1 





1 


N b + \N Y 


— M + fiVy 


K 


(0,-1,-1,0,0,0) 


3 


1 


2 

3 


1 





-1 


-N b - §AV 


M-^Ny 


bt 


(0,1,0,0,1,0) 


3 


2 


1 

6 


-1 





1 


^Vh - \N Y 


-M - \N Y 


b2 


(0,-1,0,0,-1,0) 


3 


2 


1 
6 


1 





-1 


-N b - \N Y 


M+\N y 


bt 


(0,0,0,0,1,1) 


1 


1 


-1 


-1 





1 


N b -N Y 


-M - N Y 


bs 


(0,0,0,0,-1,-1) 


1 


1 


1 


1 





-1 


-N b + N Y 


M + N y 


r+ 
c l 


(-1,-1,0,0,0,0) 


3 


1 


1 

3 





1 


-1 


N a -N b - ±N Y 


-IN Y 


c l 


(1,1,0,0,0,0) 


3 


1 


1 
3 





-1 


1 


-N a + N b + ±N Y 


iiV Y 


C 2 


(-1,0,0,0,-1,0) 


1 


2 


1 
2 





1 


-1 


N a -N b + \N Y 


\N Y 


C 2 


(1,0,0,0,1,0) 


1 


2 


1 
2 





-1 


1 


-N a + N b - \N Y 


-\N Y 


X+,Y+ 


(0,1,0,0,-1,0) 


3 


2 


5 
6 











IN Y 


IN Y 


X-,Y- 


(0,-1,0,0,1,0) 


3 


2 


5 
6 











-In y 


-§Ny 



Table 2: Different sectors and charges for the 5*0(12) model. 
3.3 Perturbative zero modes 

Given the above background, and ignoring for the time being non-perturbative effects, 
one may solve for the zero mode wavefunctions on each of the sectors of table 2. One 
obtains in this way the internal profile for each of the 4d Af — 1 chiral multiplets that 
arise from the matter curves E Q , and S c , and in particular for the 4d chiral fermions 
of the MSSM. 

Following [22] one may consider the 7-brane action derived in [5], more precisely the 
piece bilinear in fermions, and extract the equation of motion for the 7-brane fermionic 
zero modes. These equations can then be written in a Dirac-like form as 

Da$ = (3.27) 
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with 



D A = 



V 






D x 


Dy 


D z 


D x 





-D- z 


Dy 


Dy 


D- z 





-D 


D z 




D- x 






^ E 



J 



( -V2r) \ 

\ Xxy J 



(3.28) 



where the four components of ^ represent 7-brane fermionic degrees of freedom. As 
pointed out in [4], these fermionic modes pair up naturally with the 7-brane bosonic 
modes that arise from background fluctuations 

$xy = (&xy) + <Pxy = (At*) + «m (3.29) 

More precisely we have that (a^,fe) and ((p xy ,Xxy) form 4d H — 1 chiral multiplets. In 
addition, (A^, rj) form a 4d M — 1 vector multiplet that should include the gauge degrees 
of freedom of the model. One can see that these bosonic modes feel the same zero mode 
equations that their fermionic partners, and so solving (3.27) gives us the wavefunction 
for the whole M = 1 chiral multiplet. 

As we started from an 5*0(12) gauge symmetry, \l> has gauge indices in the adjoint of 
SO (12). Each covariant derivatives in D A acts non-trivially on such indices, since they 
are defined as D m = d m — i[(A m ), ■] for those coordinates m = x,y,x,y along the GUT 
4-cycle S, and as D z = —i[(Q xy ), ■} for the transverse coordinate z [22]. It is then clear 
that each sector of table 2 will see a different Dirac equation. Hence, in order to solve for 
the zero mode wavefunctions of our model, we fix the roots E p to lie within a particular 
sector of table 2 and then solve sector by sector. 

Within each sector p, eq.(3.27) is specified in terms of the following quantities: q<s>(p), 
qs(p), Qp(p) an d Qf(p)- The zero mode computation for each sector is done in detail in 
appendix A, and one can see that each of these solutions is of the form 

/ \ 



77T 



V 



m 
1 



X\D P , x\ = e-**l***-W MXty + \ y x) (3.30) 



/ 



where we have solved for the zero modes in the holomorphic gauge of eq.(3.25). Here 
are holomorphic functions of the variable \ x y + XyX, with X x , \y constants that depend 
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on the flux densities qp, qs, M xy and the mass scale m 2 , and which are different for each 
sector p. The index % runs over the different holomorphic functions that are present on 
each sector, or in other words over the families of zero modes localized in the same curve. 
Finally, recall that q<$> is a holomorphic function of the 4-cycle coordinates x, y. We have 
summarized the values of these quantities for each sector of our model in table 3, assuming 



P 



q® 



-x 



x 



-y 



x-y 



{x-y) 



X, 



A, 



A_ 



A + 

~QS\++qp 



Qs 



\--q P 



qs\+~m i 
A++<Jp-gs 



\--qp+qs 



A,-, 



A+ 

-Qs x+ - qp 



Qs 



\-+qp 



A 



+ 



A_ 



-A, 



q s \ + -m 4 

>^++qp-qs 



-A_ + 



\--qp+qs 



SU(5) rep. 



10 



10 



Table 3: Wavefunction parameters. 

for simplicity that M xy = (see appendix A for the general expressions). For the sectors 
dp , bg and c+, A + is defined as the lowest (negative) eigenvalue of the flux matrix 



-q P 

Qs 



Qp 



■ 2 \ 

im q x 

im 2 q y 



y —im 2 q x —im 2 q y 







(3.31) 



/ 



and one can check that the three lower entries of the vector in (3.30) are the corresponding 
eigenvector of this matrix. The same definition applies to A_ for the sectors a~, b q 
and c~. 5 In general ±A± satisfy a complicated cubic equation discussed in appendix A, 
which depends on the flux densities qs and qp. Since these two quantities contain the 
hypercharge flux, A± will be different for each of the subsectors , and cf. Indeed, it 
is precisely in the value of the flux densities qp and qs that the wavefunctions within the 
same SU(5) multiplet but with different hypercharge differ. 



5 Although they have a similar definition, X± have different values. Indeed, since m Q + = — m Q we 
have that — A_(a~) is the highest positive eigenvalue of m Q +. 
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Although by working in the holomorphic gauge it is not easy to see which wavefunctions 
converge locally, by going to a real gauge we can see that if we impose condition (3.16) 
the locally convergent zero modes lie in the sectors a+, 6+ and c+, see appendix A. In the 
following we will consider this choice of signs for the fluxes, and so we will concentrate in 
the wavefunctions for these sectors, that in table 4 are identified with the MSSM chiral 
multiplets. In particular, in section 5 we will compute how these zero modes are modified 
in the presence of non-perturbative effects. 



Sector 


Chiral mult. 


SU(3) x SU{2) 




q x 


% 


at 


D R 


3(3,1) 


i 

3 


-1 





at 


L 


3(1,2) 


1 
2 


-1 





bt 


U R 


3(3,1) 


2 

3 





1 


bt 


Ql 


3(3,2) 


1 
6 





1 


bt 


Er 


3(1,1) 


-1 





1 


4 


D d 


(3,1) 


1 
3 


1 


-1 


4 


H d 


(1,2) 


1 
2 


1 


-1 



Table 4: Dictionary SO(12)-MSSM 



Finally, one may solve the wavefunctions for the bulk sector (X, Y)^, which is only 
sensitive to the presence of the hypercharge flux. Although the global properties of (Fy) 
can be chosen so that no chiral matter arises from this sector [5], there will always be 
massive modes which we can be identified with the X ± , Y ± bosons of 4d SU(5) GUTs. 
As analyzed in appendix A such massive modes will have a Gaussian profile, a fact that 
can be used to suppress operators mediating proton decay [37]. 

3.4 Yukawa couplings 

Let us compute the Yukawa couplings between the chiral zero modes of this model, before 
any non-perturbative effect is taken into account. As in [22] such couplings arise from the 
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trilinear term in (2.1), which in terms of the wavefunctions above can be written as 

Yll = mJ paT [ det(^,^,^*)dvol s (3.32) 
Js 



par 



where f paT = —iTi ([E p , E a ]E T ), dvo\ s = 2u = dx A dy A da; A dy and the vectors 
are given by the three lower entries of From the last subsection and the results of 
appendix A we have that these vectors read 



m 2 



V 



Xa v 



J 



Vq 9 



V 



1 



J 



i((c r - A c 



\ 



C r 



J 



(3.33) 



where the A's and £'s are real constants defined in appendix A. The scalar wavefunctions 
X are given by 



y — m ry e (x-y)((c r X-(Xe r -(c r )y) 



x l q = m.e^-^gjix + ^y) 



(3.34) 



and following [10] the holomorphic factor can be chosen as 



ft = iiyrtv + Ca p x) 3 - 1 9j = iiy-^x + c bq yf~ j i,j = 1, 2, 3 (3.35) 

with the normalization factors 7* p , tjJ and 7 Cr to be fixed later. 

We then see that the structure of wavefunctions and Yukawas is quite similar to the 
one in the Z7 (3) toy model analyzed in [22]. One difference is the more involved sector 
structure, which as illustrated in table 4 is necessary to accommodate the MSSM chiral 
spectrum. Notice also that, due to the extra components of (F) that we have introduced, 
the holomorphic factors in the wavefunctions not only depend on the complex coordinate 
along the matter curve, but also on the transverse one. This however does not affect 
the general result of [18], in the sense that the Yukawa matrices are of rank one. Indeed, 
substituting in (3.32) shows that the integrand is the product of fcgj times an exponential 
whose argument is invariant under a diagonal U(l) rotation of x and y. Since dvol^ is 
also invariant under such rotation, the integral can be non-vanishing only when and 
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gj are constants, which happens for i — j — 3. Computing the integral yields the only 
non-zero coupling 



where f a + b + c + are the structure constants of SO (12) in the fundamental representation. 
Except for the normalization factors this coupling does not depend on the fluxes. 

4 Non-perturbative effects in local models 
4.1 The non-perturbative superpotential 

As can be seen from the previous section, the SO(12) model yields the right particle 
content but an oversimplified flavor structure. In particular, given the zero mode wave- 
functions above only one generation of down-type quarks and leptons will develop non- 
vanishing Yukawa couplings. This feature has been shown to be general for F-theory 
models where all D-type Yukawa couplings arise from a single triple intersection, and a 
similar statement holds for U-tjpe Yukawa couplings in points of E 6 enhancement [18]. 

Following [21], one may solve the above rank-one Yukawa problem by considering 
the contribution of non-perturbative effects to Yukawa couplings. Indeed, as shown in 
there, the presence of an Euclidean D3-brane instanton wrapping a 4-cycle S np of the 
three-fold base B may induce a non-perturbative correction to the tree-level 7-brane 
superpotential (2.1). 6 Such non-perturbative correction will not only depend on the 7- 
brane 4-cycle S, but also on the 4-cycle S np that the D3-instanton is wrapping, and which 
is characterized by a holomorphic divisor function h: S np = {(x, y, z) e B\h(x, y, z) = 0}. 
Indeed, by repeating the computations of [21,22] (see also appendix C) one obtains a full 
superpotential for the 7-brane on S of the form 



where the first contribution is nothing but the tree- level superpotential (2.1) while the 

second is the non-perturbative correction W np created by the non-perturbative effect. 

6 A similar effect is sourced by a gaugino condensate on a 7-branc wrapping 5 n p, see [21,49], but for 
simplicity in the following we will focus on the case of D3-brane instantons. 




(3.36) 



W — rri 



Tr($ xy F) AdxAdy + ^J2 j 6 " STr i^y F A F ) 
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Here (x,y) are local complex coordinates on the 4-cycle S (which is locally described 
as S = {z = 0}) and 9 n stand for the n th - derivative of \ogh along the holomorphic 
coordinate z transverse to S. Finally, e is a suppression factor that measures the strength 
of the non-perturbative effect. More precisely we have that 



ho = J s h the mean value of h in S and N m = (87T 2 )" 1 / s Tr(F A F) G N the total D3- 
brane charge induced on S by the presence of the magnetic field F. In addition, \i ~ m* 
is the fundamental scale for the non-perturbative effect, and A a holomorphic function of 
the closed string moduli fields which will not play any role in the following. 

Notice that the analysis of [22] did not consider the general non-perturbative super- 
potential W np above, but rather the particular case where only 9\ (denoted 9 therein) was 
non-vanishing. The reason for such Ansatz was the assumption taken in [21] that the 
two 4-cycles S and S np do not intersect, which implies that h\s = ho and so 9 vanishes. 7 
Taking 9\ as an arbitrary holomorphic function on S and neglecting all those terms with 
higher suppression on the scale m~ 2 the Ansatz of [22] follows. Now, while these are valid 
assumptions for a large class of local F-theory models (including the U(3) model of [22]) 
we will see that they need to be reconsidered for the SO (12) model of the previous section. 
As a result, the computation of zero modes and Yukawas at the non-perturbative level 
will have to be revisited to include the effects of the more general superpotential (4.1). 

In the present context, the assumption S fl S np = is well-motivated if the 4-cycle 

S localizes the MSSM degrees of freedom of the F-theory compactification. Indeed, for 

S H S'np 7^ D3-instanton zero modes charged under the MSSM gauge group may arise at 

the intersection of the two divisors. These extra zero modes would in principle invalidate 

the analysis of [21] which, similarly to [38] and [39], assumes that all instanton modes 

charged under the gauge group of interest are massive. However, it can still happen 

7 This can be understood as follows: if the intersection of two divisors Si (~1 £2 is homologically trivial, 
then the restriction of the line bundle Cs t into £2 is trivial, and vice versa. Hence, as the divisor function 
hi of Ei is a section of £e i; we can always take /ii|s 2 as a constant section of the trivial bundle £i|s 2 . 




(4.2) 



(4.3) 



with T r 



np 



Vol (S^p) + % f s C 4 the complexified Kahler modulus corresponding to S, 
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that S PI S np 7^ if S is wrapped by a 7-brane that does not localize the MSSM gauge 
group. In fact, from the results of [40] we know that [S fl S np ] 7^ must be true for 
at least some 7-brane 4-cycle S, since otherwise a D3-instanton wrapping S np would not 
have the right structure of zero modes to generate a superpotential. This result is in 
perfect agreement with the intuition of D3-instanton generated superpotentials in type 
IIB orientifold compactifications. Indeed, if our F-theory compactification has a type IIB 
limit, then by taking it the 4-cycle 5* such that S fl S np 7^ will become the 4-cycle of an 
07-plane, while the D3-instanton zero modes localized at S fl S np will correspond to the 
universal neutral zero modes of an 0(1) D3-brane instanton [41]. 




Figure 1: Two possible scenarios for non-perturbative corrections to Yukawa couplings. In figure 
i) the D3-brane instanton intersects a 7-brane that in turn intersects the 4-cycle Squt = Smssm- 
In figure it) the D3- instanton intersects only 7-brane divisors fully disconnected from Squt- 

To summarize, even if the D3-instanton generating the non-perturbative superpoten- 
tial W np within (4.1) does not intersect the 4-cycle Smssm that localizes the MSSM degrees 
of freedom, it must intersect some other 7-branes and this, in a potential type IIB limit, 
should correspond to a transverse intersection of the D3-instanton with an 07-plane. 
Given this, one may consider two different scenarios: 

i) The D3-instanton intersects a 4-cycle S' which in turn intersects the 4-cycle Smssm 
in a matter curve S. That is, [Sd3 H S'] 7^ 7^ [S' fl Smssm] but [Sd3 H Smssm] = 0. 

ii) The D3-instanton intersects a 4-cycle S' that does not intersect Smssm at all. That 
is, [S D3 nS']^0 and \S' n S MS sm] = = [S D3 n Smssm]- 
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Let us consider both scenarios in the context of the 5*0(1 2) model of the previous 
section, taking advantage that it can also be understood as a local type IIB orientifold 
model. Indeed, if we interpret the 5*0(12) model in terms of D7-branes at angles, we 
obtain a configuration of the form 



Brane 


Stack 


Divisor 


07-plane 




S : z = 


1 D7-brane 


a 


S a : z = a(\y - x) 


1 D 7- brane 


a* 


S a * : z = -a{\y-x) 


5 D7-branes 


P 


Sp: z = \ay 


5 D7-branes 


(3* 


Sp* : z = -\oy 



where we have ignored the world volume flux F on each stack. Here we have defined 
the adimensional quantity a = (m/m st ) 2 , with m~ t 2 = 2ira' being the string scale in 
the orientifold limit, that determines the angle of intersection between D7-branes. The 
D7-branes (a, a*) are mapped to each other by the orientifold action, same for ((3,(3*), 
and so we only get a unitary gauge group from each pair. Namely, we obtain a gauge 
group £7(5) x U(l) a ~ SU(5) x U(l) a x U(l)p, as expected. Finally, we can also match 
the matter spectrum obtained in the previous section by realizing that the sectors a, b, c 
arise from the matter curves S a = S a fl Sp, S fe = Sp fl Sp* and S c = S a fl Sp, and then 
including the effect of the world volume flux. 

In this type IIB version of the 50(12) model it is clear how to realize the scenario 
i) above. First we should require that the 4-cycle 5 np wrapped by the D3-instanton 
intersects the 07-plane 4-cycle 5, and then that it does not intersect S a , Sp or their 
orientifold images. Notice that by intersecting the 4-cycle 5 we do not generate any extra 
D3-instanton zero mode. On the contrary, this nontrivial intersection is necessary for 5 np 
to host a 0(1) D3-instanton and so to get rid of unwanted zero modes. 

Hence, in this scenario we have that the D3-instanton intersects the 07-plane in a 
non-trivial 2-cycle Se3 = 5 fl 5 np , which by assumption is far away from the Yukawa 
point 5 fl S a fl Sp at x = y = z = 0. By the discussion above, this means that the 
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pull back h\s of the divisor function h will not be a constant, nor will be the function 
6 oc [log(/i//io)] 2 =o that enters into the superpotential (4.1) if 5 refers to the 4-cycle of 
the 07-plane. Now the important point is that, for the SO (12) model under discussion, 
the 4-cycle S in (4.1) indeed corresponds to the embedding of the 07-plane, rather than 
to Smssm- 8 As a result, for the 5*0(12) model and within the scenario i) above, 9 (x,y) 
will be a non-constant holomorphic function on S. 

On the other hand, if we consider the scenario ii) for the SO (12) model we have that 
S H 5 np = 0, and so in this case 9 = just like in the Ansatz of [22]. In fact, as we 
will see momentarily, it happens for the 50(12) model that the contribution to (4.1) 
coming from d\ also vanishes, and so for this scenario the first non- vanishing contribution 
to the non-perturbative Yukawa couplings would come from 2 , which is already rather 
suppressed in the scale m. Thus, for the purpose of generating realistic Yukawa couplings, 
the previous scenario seems more promising. 

We then see that the two scenarios i) and ii), which have a clear geometrical meaning 
in terms of a global model (cf. fig. 1), can be characterized in terms of 6 in our local 
description. Such function should be taken constant for scenario ii), while it will be given 
by a certain holomorphic function for scenario i). Finally, this result will also hold away 
from the type IIB orientifold description of the 50(12) model, the only difference being 
that we should replace the 07-plane by a more complicated system of (p, q) 7-branes. 

The discussion above shows how the geometry of our F-theory compactification can 
constrain the superpotential (4.1), and one may find further geometrical constraints on 
these non-perturbative corrections. Indeed, let us go again to the type IIB orientifold 
description of our local 50(12) model. There the non-perturbative piece of the super- 
potential (4.1) is generated by a 0(1) D3-instanton, whose embedding 5 np = {h = 0} 
determines the holomorphic functions 6 n via the divisor function h. Now, the 4-cycle 5 np 
will host a 0(1) D3-instanton only if it is invariant under the orientifold action, which 
close to our Yukawa point is given by z i— > —z. This means that in the vicinity of this 
point we have to impose h(z) = h(—z), and from the definition (4.3) it follows that 9 n = 

8 Note that for ($) = we recover an SO (12) gauge group localized at the 4-cycle S, namely a stack 
of 12 D7-branes on top of the 07-plane locus, while for ($) ^ the gauge group is broken because these 
D7-branes are wrapping different 4-cycles, one of them being Smssm- Hence Smssm S. 
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for n odd. In particular the function 9\, in which the Ansatz used in [22] was based, does 
vanish for the SO (12) model, and so in principle the whole analysis needs to be revisited. 

Again, this kind of constraint is not only true in the type IIB orientifold limit, and 
it holds for general local F-theory models. In fact, one can see directly from the super- 
potential (4.1) that the terms that correspond to 9 n with n odd vanish automatically for 
SO(N) groups, in agreement with our type IIB intuition. Indeed, such terms are propor- 
tional to symmetrized traces of the form STr(tit2 • • - i n +2), with tj generators of SO(N) 
in the vector representation. As these tj are antisymmetric matrices, the symmetrized 
trace will vanish identically for any n which is odd. In particular, the symmetric tensor 
d abc = STr(t a t 6 t c ) which corresponds to the correction of 9 1 and played a key role in [22] 
will also vanish identically. As already pointed out in [22], the fact that d abc = is not 
only true for SO(N), but also for other groups of interest in F-theory GUT model building 
like E 6 , E 7 and E 8 . Hence, an analysis of non-perturbatively generated Yukawas that not 
only involves 9\ ^ but the more general superpotential (4.1) is in order. 

4.2 The corrected equations of motion 

As a consequence of the new superpotential (4.1), the equations of motion for the 7-brane 
fields $ and A are modified, affecting a local model in two different ways. On the one hand 
the background values (<&) and (A) will be shifted to a new value. On the other hand the 
fluctuation fields (a m , ipm) and (<p xy , Xxy) will satisfy different zero mode equations, and so 
the wavefunctions obtained in the previous section will receive a correction that accounts 
for the non-perturbative effect. In the following we will discuss the new equations that 
the 7-brane fields need to satisfy, solving for the background to first order in the non- 
perturbative parameter e. We will also analyze the new zero mode equations for the 
bosonic fluctuations (a™, <Pxy), which by supersymmetry satisfy the same equations as the 
fermionic modes (ipm,Xx y )- The next section will be devoted to apply this latter set of 
equations in order to obtain the corrected zero mode wavefunctions for the SO (12) model 
of section 3. 
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The equations of motion that follow from the superpotential (4.1) are 
DA y + eJ2 {D y [S(0n® n xy Frn X )} + D x [S(9 n $ n xy F yfh )} + D^[S(9 n ^ y F xy )]) = (4.4a) 

n 

F A dx A dy + e - J2 n0 n S(® r xy lF A F) = (4.4b) 

n 

where S stands for the symmetrization of all the elements within the brackets, and D m 
are the covariant derivatives along the GUT 4-cycle defined below eq.(3.28). Notice that 
if we take all 9 n = except for 9\ we recover eqs.(3.40) of [22], as expected. 
In addition, the D-term (2.2) gives the additional equation 

wAF+J[$,f] = (4.5) 

where in the local coordinate system that we are working with the Kahler form is taken 
to be 

uu — — {dx A dx + dy A dy) . (4.6) 

Notice that, unlike the F-term equations (4.4), the D-term equations do not receive any 
correction due to the non-perturbative effects. This may seem surprising, but it can be 
derived directly from the results of [42], as we discuss in appendix C. 

While the system of differential equations (4.4) and (4.5) is quite difficult to solve, one 
may follow [22] and perform a perturbative expansion in the small parameter e 

**y = *S + € Sg) + e 2 Sg> + • • • A m = A^+eA^ + ^A^ + ... (4.7) 

and then solve for the above equations order by order in e. Indeed, to zeroth order in this 
expansion (4.4) reduce to 

(An<M (0) = *& = (4-8) 

which are indeed the 7-brane F-term equations in the absence of any non-perturbative 
effect. To first order we obtain 

(D r -A y p = - (d y 9 )F^ x - (d x 6o)F$ 

" E n {Sl^Dyie^F^} + S^D^e^F^ (4.9a) 

n=l 

F« AdxAdy = -^n^S^'FAFf ' (4.9b) 

n 
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where we have used the Bianchi identity D[ m F np ] = and we have defined 

= D<®*V-i[A$ i \*W] (4-10) 

Let us now provide a solution to these 0(e) equations at the level of the background, 
generalizing the results in eqs.(3.45) of [22], where they were solved for the particular case 
where 6 n = S n i9. 9 This more general solution is given by 

(A^> {1) = ~$>*» S [(^) n_1 ((^)(^) - (A y )(F xfh ))] {0) (4.12) 

n=l 

which indeed solves (4.9b) at the level of the background. Note that this result is true 
even if the matrices in (4.12) do not commute. Similarly, we have that 

(<s> xy )W = ((A y )^d x e -(A x )^d y e ) (4.13) 

+ \ j2 n ( S l(A y )(® xy ) n -\d x + D x ){6 n (<$> xy ))\ 

n=l 

- S[(A x )^ xy ) n -\d y + D y )(9 n ^ xy ))}) {0} 

satisfies eq.(4.9a) at the level of the background. Indeed, by applying = d m to the 
rhs of (4.13) and using that 6 n and (Q xy )^ are holomorphic we obtain the rhs of (4.9a) 
plus the term 

^£ n£ # (S [(^) n " 1 (^)[(^), <***>]] ~ S [(^) n - 1 (A>[(^->, (^)]]) (0) 

n 

= \H n()l n P , (^)] (0) = i [(A*)®, (^) {0) ] (4.14) 

n 

as needed from eq.(4.9a). Here i,j — x, y and 9* y = —Q y * = —0 n . 

9 As pointed out in [22] these solutions for the background are directly related to a Seiberg-Witten 
map [43] acting on the 7-brane fields A and $. Hence, a simple way to find the solutions for (4.9) is to 
generalize the SW map in [22], eqs.(4.18) to the present case. A natural generalization is given by 

A m = Am +A ih =A ih - e -Y j n 0% S [^MdjA^ + F jm )] + 0(e 2 ) (4.11a) 

n 

K - + *-.=*-y-|E S \- Ai ^ + D iWn^)] + 0(e 2 ) (4.11b) 

n 

where i, j = x, y and 6* v = — 0^ = —6 n - At the level of the background and in the holomorphic gauge 
(A m )W = 0, we can identify (A m ) = -e(A m )^ and {$ xy ) = -e{$ xy ) {1 \ obtaining eqs.(4.12) and (4.13). 
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The next step is to derive the zero mode equations for the bosonic fluctuations 
(<Pxy, «m) defined by 

<5>xy = (®xy)+(Pxy A m = (Afn) + dm (4.15) 

and which also follow an expansion of the form (4.7). Expanding the O(e ) equation (4.8) 
to linear order in fluctuations we obtain the equations of motion for ((f x ° y \ a^)- O ne can 
check that together with the zeroth order D-term equation, they amount to the system 
of equations (3.27), where now ^ = (0,a^\af\(p x ° y ). Hence at zeroth order in e the 
wavefunctions for the bosonic fluctuations match the fermionic wavefunctions in absence 
of non-perturbative effects, as expected from supersymmetry. For the 5*0(12) model, the 
zeroth order wavefunctions for each sector are then given by (3.33). 

Let us now consider the equations for (<p£y , a^), that arise from expanding (4.4) to 
first order in e. Just like in [22], in order to write down the zero mode equation to order 
O(e) we need to take into account the corrections to the background (4.12) and (4.13). 
Indeed, let us first consider the particular case where 9o ^ but 9 n = for all n > 0. 
Then, at the linear order in fluctuations and in the holomorphic gauge, eqs.(4.4) read 

D x a y — D y a x = (4.16a) 
DrnVxy + i[($xy), Ofh] = * [(d y 6o)D x a m - (d x 9o)D y a m } (4.16b) 

In addition, the corrected background solutions reduce to 

= + O(e 2 ) (4.17) 
($^> = ($ xy )(V + e((A y )^d x 9 -(A x )^d y 9 )+O(e 2 ) (4.18) 

and so, by plugging the latter into (4.16) we obtain 

d x ay-d y a x = + O(e 2 ) (4.19) 
drn<p xy + i[(<S> xy )(°\arn} = e [(9^)^ - (^o)^] + 0(e 2 ) (4.20) 

which could have also been obtained by linearizing (4.8) and (4.9) in fluctuations. Notice 
that these F-term equations are independent of the worldvolume flux, and that the cubic 
coupling that arises from (4.1) for only 9 ^ is also flux-independent. Hence, following 
the arguments of [18] we have that the holomorphic piece of the Yukawa couplings should 
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not depend on these fluxes. Moreover, it is easy to see that shifting 9 by a constant will 
not modify nor the wavefunctions nor the Yukawa couplings. This can be traced back to 
the fact that a constant 9 will only add a constant piece to the superpotential (4.1). A 
detailed analysis of the F-term zero mode equations for this case as well as a computation 
of the holomorphic Yukawas in terms of a residue formula is given in appendix D. 

Similarly, the zero mode equations for (ifxy,a^) with general 9 n can be obtained by 
linearizing eqs.(4.9) in fluctuations. Let us do so for the case where only 9 and 9 2 are 
non-vanishing, which will be the case analyzed in the next section. From (4.9a) we have 

^g, ) +t[($xy> (0) ,a2 ) ] = ( d y e )d x a^ - (d x 9 )d y a^ (4.21) 

+ 2S[(<f> xy ) (d y (9 2 ($ xy ))D x a m - d x (9 2 (<f> xy ))D y a m )]W 

+ 2% l($ xy ) (d y (9 2 (<f> xy ))(A x ) - d x (9 2 ($ xy ))(A y )) , a m f ] 

+ 2S[($ xy ) (D y (9 2 <p xy )(F xm ) - D x (9 2 <p xy )(F ym ))}^ 

+ i9 2 [{&xy) ((F X m)(A y ) — (F ym )(A x )) 

+ 2S[<p xy (d y (9 2 (<f> xy ))(F xm ) - d x (9 2 (<S> xy ))(F m ))]W 
where we have used the assumption that (Am) and (& xy ) commute. From (4.9b) we obtain 

d x af - d y a x 1] = 29 2 S [cp xy ((F xx )(F yy ) - (F^)(F„))] (0) (4.22) 

+ 29 2 S [($ xy ) ((F xx )D y a y - (F yx )D x a y + (F yy )D x a x - (F xy )D y a x )f ] 

+ %9 2 [($ xy )((F xx )(A y ) — (F yx ) (A x )) , a y ]^ 

+ i9 2 [(& xy )((F yy ){A x ) — (F xy ) (A y )) , a x }^ 

which for non-vanishing fluxes are rather involved equations. Nevertheless, by means 
of the Seiberg-Witten map of footnote 9 one can simultaneously get rid of all the flux 
dependence within these two F-term equations and in the cubic coupling that arises from 
(4.1) with 9 2 7^ 0, see appendix D for details. Hence, according to the results of [18] 
we have that the holomorphic piece of the Yukawa couplings should not depend on the 
worldvolume fluxes, a fact that we will use in the next two sections. 

Finally, the fluctuations must satisfy the D-term equation that can be derived from 
expanding (4.5) to linear order in fluctuations, namely 

u>Ad {A) a-±[(Q), ( p] = (4.23) 
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where a = a x dx + dydy and tp = <p xy dx A dy. As mentioned above, at zeroth order in e we 
obtain part of the set of equations (3.27) for (ip x °y,a!^). At first order in e we obtain 



D x af + D v a$) -i[fe) (0) ,^] = < [(^) (1) ,4° } J + i [<AA« afj +< [<^> (1) >¥>S] 

(4.24) 

where (A m )W and (®xy)^ are the Hermitian conjugates of (4.12) and (4.13), respectively. 



5 Non-perturbative zero modes at the SO (12) point 

Let us now apply the non-perturbative scheme of the previous section to the Yukawa point 
of 50(12) enhancement. More precisely, we would like to consider the superpotential (4.1) 
for the local £0(12) model described in section 3, and see how this new superpotential 
affects the wavefunction profile for the matter fields and the Yukawa couplings. The aim 
of this section is to solve for the zero mode wavefunctions in the presence of the non- 
perturbative piece of the superpotential, leaving the computation of Yukawa couplings for 
the next section. 

Before attempting to solve for such wavefunctions, it is useful to recall some results 
obtained in the previous section, which impose some constraints on the superpotential 
(4.1) for the model at hand. Indeed, because the 7-brane fields $ and A take values in the 
algebra of 50(12), all those non-perturbative contributions in (4.1) with n = odd vanish 
identically. Hence, we should consider those terms that involve the holomorphic functions 
9 , 9 2 , ... etc. 

As also discussed in the previous section, the new zero mode wavefunctions can be 
written as a perturbative expansion in the small parameter e that measures the strength 
of the non-perturbative effect. More precisely we have that 

->■ -Ho) /-Hi) -Hi) \ 

^ = i/+^Kk + i P k + ---)+^ 2 ) (5-1) 

—Ho) 

where ip are the tree-level wavefunctions (3.33)-(3.35) for the sector p = ap,6+,c+, 
-Hi) 

■0p,0 o is the 0(e) correction to this sector when all #2n vanish except 6q, and the same for 
-Hi) 

ip 2 . If we consider both 6 and 6 2 present at the same time, the first order correction to 
the zeroth order wavefunctions must satisfy the F and D-term equations (4.21), (4.22) and 
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(4.24). Notice that 9 and 9 2 appear linearly in the rhs of these equations, which justifies 

-Xi) -Xi) 

that we can write the first order correction as the sum vp p 6(j + ip p02 . This statement 
generalizes to the case where more 9 n are present in the non-perturbative piece of the 
superpotential. However, as we will see in the next section, for the purpose of computing 
Yukawa couplings considering only 6> and 9 2 is enough. 

-Xi) 

In the remaining of this section we will solve for the first order corrections p % n to 

the wavefunctions (3.33)-(3.35), first for the case where n = and then for n = 2. As it is 

clear from eqs.(4.21)-(4.24) the corrections for the case n = 2 are more difficult to obtain 

but, as mentioned before, one can do a field redefinition that removes the worldvolume 

-Xi) 

flux dependence. Hence, in section 5.3 we will determine the corrections ip Pj g 2 with fluxes 
turned off. Still, the computation for the n = 2 case is slightly more technical and the 
reader not interested in such details may focus on the simpler n = case in order to get 
the idea of the computation, and then proceed to the computation of Yukawa couplings 
in the next section. 



5.1 Zero modes for n = 



Let us then consider the wavefunction corrections for the case where 9q ^ but 9 n = 
for all n > 0. Notice that in [22] 9 was also present but assumed to be constant, and 
shown that Yukawa couplings are independent of it. As discussed in the previous section 
we may now relax this condition and take 9q to be non-constant and holomorphic on x, y. 
For simplicity let us take it to be 



9 = im 2 (6»oo + x 9 0x + y 9 0y ) 



(5.2) 



where 9 00 , 9 0x , 9 0y are complex constants, and the factor im 2 has been added for later 
convenience. We will now see that the corrected wavefunctions do depend on 9 0x and 9 0y , 
and in the next section that they also enter into the corrected Yukawa couplings. 

The corrected equations of motion for the fluctuations were derived in section 4, and 
they can be conveniently rewritten by using the notation 



/ „_ \ 



(p) 



''px 



u py 

y tPpxy J 



(5.3) 
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Recall that because of supersymmetry the bosonic fluctuations (a P m, f P x y ) pair up with 
fermionic fluctuations (ip P m, Xp) analyzed in section 3.3, and so in the absence of non- 
perturbative effects the components of ip match those of the vectors in eq. (3.33). Using 
the above notation we can express the corrected F-term equations (4.19) and (4.20) as 

dscipy - dydpx = (5.4) 
dfnVpxy + im 2 q^(p)a P rn = iem 2 [9 0y d x aprn - Ox d y a P rn] + 0(e 2 ) (5.5) 

which, in the holomorphic gauge that we are considering, do not depend on the world- 
volume flux densities M xy , qp and q s . By the results of [18,20] this implies that the 
holomorphic Yukawa couplings cannot depend on these quantities either. We will verify 
this explicitly by means of a residue computation performed in appendix D. 
On the other hand, the D-term (4.24) translates into 

{d x + x[q P (p) - M xy q F (p)\ - yq s (p)} a px 

+ {d y ~ y[qp(p) + M xy q F (p)} - xq s (p)} a p y - im 2 q^(p)ip pxy (5.6) 
= iem 2 9 0x {y[M xy q F (p) + qp{p)\ + xq s (p)} p pxy 

- iem 2 9 0y {x[M xy q F (p) - q P (p)} + yqs(p)} ¥ P xy 

where the specific values of the charges q®, q F , qp, and qs, for each sector p are given in 
tables 1 and 2. Notice that as usual the D-term equation depends on the flux densities, 
and in particular on the hypercharge fluxes contained in qp and qs- Hence, just like in [22], 
the holomorphic Yukawa couplings will not depend on the hypercharge but the physical 
Yukawas will, as we show in the next section. 

As already mentioned, the zero modes to zeroth order in e are given by eq. (3.33). To 
find the corrected zero modes the strategy is to start with an Ansatz motivated by the 
zeroth order solutions and then proceed perturbatively in e. Notice that the zeroth order 
solutions consist of a fixed vector lt p given by (A. 54) multiplying a scalar wavefunction 
Xp = Vpxy, that has a simple dependence on the complex variables X x y+X y x and X x x—X y y. 
We find that the first order solutions are also of this form, but now with a corrected scalar 
wavefunction, namely 

<P P *y = A + VP™ + 0{e 2 ) (5.7) 
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where tpf xy are given by the scalar wavefunctions Xp m ec f (3.34). In the following we 
report the results for the correction <f$ y in the different sectors, dropping the subscripts 
xy to avoid cluttering the equations. While we will keep the worldvolume flux dependence 
in the zero mode equations, for simplicity we will set M xy = 0. Our solutions are however 
easily generalized for non- vanishing M xy . 

Sector a + 

Let us first define the complex variables 

U a = X-CaV Va = V + CaX (5.8) 

that come from a rescaling of X x y + XyX and \ x x — \y-y for this sector. To simplify notation 
we have suppressed the subindex p = 1,2, for each sector p = ai", that labels elements of 
the 5 representation with different hypercharge, and therefore with different values of A a 
and ( a - Notice that the variables u a and v a are different for each of these subsectors. 

The correction to the wavefunctions ip l a+ can be conveniently written in terms of 
(u a ,v a ). Concretely, 

^ = (Mv a ) + T a+ ) (5.9) 

where g$(a + ) = — x for these sectors, Ai is a holomorphic function of v a and 



T a+ = ^Oy^l U lfi( V a) + - 0O*)//O>a) + y^ + W fi(v a ) (5.10) 

The last two terms are in fact only necessary to fulfill the D-term equation (5.6). In 
particular, we obtain that the coefficients ct\ and a 2 must take the following values 

4 

TfL — — 

«i = -^{0oA<ls-(aq a p)+0oy(q a p + (aq a s)} 
A a 

4 

TfL — — 

«2 = - t {Mp-W} (5.H) 

To sum up, taking into account the zeroth order solution in eq.(3.34) the final result 

for (f l a+ can be cast as 

rt + = m*e-«* A ^ (f t (v a ) + eA t (v a ) + eT^ + ) + 0(e 2 ) (5.12) 

Naively the functions Ai remain unfixed by the above equations of motion. This is because 
we could think of them as a 0(e) correction to the holomorphic functions j\, which are 
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also not fixed. However, given the choice (3.35) of the Ai are fixed as follows. Notice 
that the F-term equation (5.4) implies a pin = <9m£ p , where £ p is a regular function [18]. 
As shown in appendix D, this function can be found by integrating (5.5). Imposing that 
£ p is regular at the loci q$(p) = where the zero modes are localized implies nontrivial 
constraints for the wavefunctions. In particular, for the a + sector requiring to be free 
of poles at x = fixes the A i: which read 

A 2 = A 3 = ; A 1 = 1 1 a a ; a = ml( a (( a 9 0y - 29 0x ) (5.13) 

Interestingly, this form of Ai guarantees that the Yukawa couplings computed via overlap 
of zero modes will be flux independent up to normalization factors, as we comment on 
the next section. 

Sector b + 

The sectors b + and a + are quite similar, so let us first define the variables 

u b = y-( b x Vb = x + ( b y (5.14) 

that again are different for each subsector p = &+, q = 1, 2, 3. As in (5.9) we obtain 

X* 1 ) = m,e«* A ^ (Bj(v b ) + V b ) (5.15) 
where now q<s>(b + ) = y, Bj is a holomorphic function of v b and 

= ^0ox^l u b9j(vb) + hu b ((b6ox - 6oy)g'j(vb) + 

Again, the terms proportional to fii, @2 are required to satisfy the D-term equation. These 
coefficients can be obtained from (5.11) by replacing 9 0x <H> 9 0y , qp — > —q b P , q% — > q b s , and 
consistently A a ->• X b , ( a ->■ ( b . 

Including the zeroth order solution from eq.(3.34) yields the full result 

rt + = m^* x ^ ( 9j (v b ) + eBjivb) + eT{) + 0{e 2 ) (5.17) 

The functions Bi are determined demanding that be free of singularities as explained 
in appendix D. In this way we find 

B 2 = B 3 = ; B l = 1 lb ; b = m%{CbO Gx - 29 0y ) (5.18) 
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\y 2 + $2yv b 



9j(v b ) (5.16) 



Sector c 



In this case it is helpful to introduce the variables 

u c = x-r c y v c = y + r c x (5.19) 

with r c = (A c — Cc)/Cc- The quantities A c and ( c , defined in appendix A, actually depend 
on the subsector c+, r — 1, 2. The correction to the wavefunction is found to be 

<p { $ = m* lc e q ^ (C{v c ) + T«) (5.20) 

where now g$ = (x — y), C is a function of t> c , and 

T« = -\£ul% x + 9 0y ) + 6 -l( x -yf + 5 2 {x - y)v c (5.21) 

The constants 5\ and <5 2 are given by 

4 

61 = ^J]^r^2 {^x (Qs ~ r c q c P ) + 0y (q c P + T c q c s )} 

4 

777 _ _ 

62 = ai - Te)2 {gox (gp + <f s ) + goy (gp ~ qs)} (5-22) 

The holomorphic terms in , which depend on 6q x and 9q v through 8\ and 82, are needed 
to satisfy the corrected D-term equation. In appendix D we show that C = 0. 

5.2 Normalization and mixings of corrected zero modes 

Given the corrected zero mode wavefunction one must compute their normalization fac- 
tors, since it is through these factors that the physical Yukawa couplings depend on world- 
volume fluxes. The computation of normalization factors for perturbative zero modes is 
carried out in appendix A, where the following norms are computed explicitly 

— Vreal — Vrcal „ f — Vrcal — Vtreal . 

K; J = {1> pi \1> pj ) = ml Tr(^ • i> pj )dvol s = 2||^ p || 2 X^ (5.23) 

Js 

with if p defined in (A. 54), and the scalar wavefunction metrics 

X'J = ^* (Xp real > 4 reaI ) = ml [ (x" ca T4 rcal ^ol 5 (5.24) 

Js 

are calculated by extending the local patch to C 2 . Here the superscript 'real' stands for 
the zero mode wavefunction expressed in a real gauge rather than in the holomorphic 
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gauge that we have used in the previous sections. Following [11] one may switch from 
wavefunctions in the holomorphic to the real gauge by multiplying them by an appropriate 
sector dependent prefactor. For instance, in the sector a + we have 



i real 



where x\+ is the zero mode computed in the holomorphic gauge. As we are now deal- 
ing with non-perturbative zero modes, we should take x\+ to be the corrected scalar 
wavefunction ip l a+ given in (5.12). 

For the perturbative zero modes analyzed in appendix A, X l J = for i ^ j, since 
the integrand needs to be invariant under the U(\) diagonal rotation (x,y) — > e ta (x,y). 
However, this no longer needs to be true at order e. Indeed, let us consider the sector a + 
and the corrected zero modes in this sector when only 9q is present. Due to the correction 
given in (5.10), one can in principle have non-diagonal metrics. In fact, to order e 
one finds that only and its conjugate are different from zero. 

Substituting the corrected zero modes and computing the Gaussian integrals we find 

K + = 2 -^0a\\ 2 Y;il^ + O(e*) (5.26) 
where A a = — (2A a + q a P (l + \\~Va\\ is given by (A. 62) and the matrix x a is 



( 2^ e[ml(2r a e 0x + r 2 a e 0y )}\ 

0^0 

1p 

\e[ml(2r a 9 0x + r 2 a 9 0y )] 1 J 

-.a 



(5.27) 



" - -| (5.28) 
The quantity r a is the quotient between the off-diagonal and diagonal worldvolume fluxes 
felt by this sector. Hence, when off-diagonal fluxes are turned on is non-zero. Finally, 
notice that the diagonal terms K™ + do not get corrections to order e. 

As the metric for the zero-modes is non-trivial, the Yukawas computed from them do 
not yet correspond to the physical couplings. From the 4d effective theory viewpoint, to 
obtain physical Yukawas one performs a field redefinition such that the 4d chiral fields have 
canonical kinetic terms. The higher dimensional counterpart of this field redefinition is to 
take a linear combination of zero modes such that the matrix K a + becomes the identity. 
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In the absence of non-perturbative effects e = and such redefinitions are rather easy to 
perform, since K a + is diagonal and one only needs to choose the normalization factors as 
in (A. 65) to have K a + = 1 3 . The same applies to order (9(e) whenever = 0. 

In general we will have 4d fields $™ whose metric K p has off-diagonal terms, and so a 
wavefunction rescaling is not enough to have canonically normalized fields. However, as 
the field metric K p must be Hermitian and definite positive, we can always write it as 



K = P^P 



(5.29) 



and describe the fields with canonical kinetic terms and the physical Yukawa couplings as 



^phys — (-fp)m^ ^ijk ~ i^p )i (Pp' )j (Pp" )fe^mnp 



(5.30) 



where Y mnp stand for the initial set of Yukawa couplings. Finally, the set of internal 
wavefunctions ip m associated to these fields will transform under this change of basis as 



phys 



(p-^r^. 



(5.31) 



In the case of K a + one can easily find a matrix P a + such that (5.29) is satisfied 



Pa+ 



V2-Kml\\~^ c 



( V2^ 



\ 



(g°)V2 



/ 



la 



ll 



\ 



\ 



ll / 



+ C(e 2 ) (5.32) 



a 

Ha = -?=(2r a 9 0x + r 2 a 9 0y ) 
Upon choosing the normalization factors 7* as in (A. 65) this matrix simplifies to 



(5.33) 



P a+ ^% 



fl 



\ 



Wa ^ 



/ 



+ (9(e 2 



(5.34) 



Hence, in order to have canonically normalized fields we not only need to choose appro- 
priate normalization factors, but also perform a rotation among the families of this sector. 
One can express such rotation in terms of the holomorphic representatives fi that describe 
each of our families as 



fl — fl /2 — }2 fz — J3 — t^afl 
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(5.35) 



fi being the holomorphic representatives that describe canonically normalized 4d fields. 

The analysis of the metrics in the b + sector proceeds along the same lines. Again con- 
sidering the case where only 9 ^ 0, the non-zero off-diagonal entries of the mixing matrix 
amount to X^l and its conjugate X™. Similarly to the sector a + there are no (9(e) cor- 
rections to the diagonal terms X£+, and obtaining canonically normalized fields amounts 
to appropriately choosing the normalization factors 7^ and performing a redefinition of 
chiral representatives. More precisely, we must choose 7^ as in (A. 66) and perform the 
redefinition 

9i = 9i h = 92 h = 93- efib9i (5-36) 

where now 

fjL b = ^(2r b 6 0y + r 2 b 6 0x ) r b = % (5.37) 
V 2, q P 

To summarize, in the presence of non-diagonal fluxes of the form (3.12) non-trivial 
corrections appear at first order in e for the wavefunction metrics. The corrections appear 
in the off-diagonal metrics , i 7^ j that vanish at zeroth order in e. One can perform 
a change of basis of the families of zero modes in order to set the off-diagonal entries 
of K p to vanish, and choose the appropriate wavefunction factors 7* to have canonically 
normalized 4d kinetic terms. Notice that such factors are the ones found in appendix 
A, namely the normalization factors at the perturbative level, and that the only effect 
of (9(e) corrections to the field metrics is encoded in the redefinitions (5.35) and (5.37). 
These redefinitions have however no effect for the physical Yukawa couplings, at least at 
the level of approximation that we are working. 

Indeed, the above redefinitions are only nontrivial in the cases fa — > f 3 and g 3 g 3 , 
and amount to say that in this new basis the third families of both sectors a + and b + have 
a contamination of 0(e) from the first family. In principle this contamination will modify 
the Yukawa couplings Y 3 j and Y i3 . However, as the Yukawa couplings involving the first 
family are already of order e, the modification will be 0(e 2 ). The same result is obtained 
by directly applying (5.30) to the Yukawa couplings computed in the next section. 

We then find that, although non-perturbative effects modify the metrics for the 4d 
matter fields, this modification can be neglected at the level of approximation that we 
are working, at least for the purpose of computing fermion masses. The whole effect 
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of wavefunction normalization is already captured by tree-level wavef unctions, and more 
precisely by the computations carried out in appendix A. This fact is quite relevant in 
the present scheme, because holomorphic Yukawas do not depend on worldvolume fluxes. 
Hence, the only place where hypercharge flux will enter into the expression for the physical 
Yukawa couplings will be via the normalization factors of perturbative zero modes. 

While the above results have been obtained for the superpotential (4.1) with only 
6 7^ 0, they hold more generally. Indeed, one can check that a similar result is obtained 
for the U(3) model of [22], that corresponds to the case n — 1, and one expects the 
same for the case n — 2. That is, one expects 0(e) corrections to the off-diagonal metric 
elements, which can nevertheless be absorbed into a redefinition of the third family and so 
do not affect the expressions for the physical Yukawa couplings. In particular, the physical 
Yukawas are expected to depend only on the O(e ) normalization factors 7* computed in 
appendix A, where the full hypercharge flux dependence will be contained. 

5.3 Zero modes for n = 2 

-Hi) 

Let us now turn to compute the corrections ^ Pi e 2 defined in (5.1). To this end we 

set 6q = 0, while for simplicity we take 9 2 ^ to be constant. Note that the F-term 

equations given in (4.21) and (4.22) are rather unwieldy when fluxes are present, but 

simplify considerably when fluxes are turned off. On the other hand, just like in [22] and 

in the case n = above, the holomorphic Yukawa couplings for the case n = 2 are expected 

to be independent of worldvolume fluxes. We show that this is the case in appendix D, by 

means of the Seiberg-Witten map of footnote 9. Indeed, it is easy to see that the SW map 

(4.11) defines certain variables and ip xy for the 7-brane fluctuations. When we express 

the F-term zero mode equations and the superpotential trilinear couplings in terms of 

such hatted variables, we find that they are flux-independent. 10 We can then compute 

the holomorphic Yukawa couplings via a residue formula as in [18], showing explicitly 

their flux independence. 

Since we know that holomorphic Yukawas are flux independent, in order to compute 

them we may solve for such wavefunctions by considering the zero mode equations in the 

10 More precisely, the F-term equations and trilinear couplings for <p xy ) are those for (a™, (p xy ) 
after setting (F) = 0. 
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absence of worldvolume fluxes. Of course, we are ultimately interested in computing the 
physical Yukawa couplings, which do depend on fluxes. However, as we have argued in 
the last subsection, this flux dependence is expected to arise only via the normalization 
factors for the tree-level zero modes. Hence, one can simplify the computation of Yukawas 
by switching off worldvolume fluxes when computing corrected wavefunctions and their 
overlap integrals. The flux dependence is taken into account by including the normaliza- 
tion factors 7* calculated in appendix A. This last step will be carried in the next section, 
where the Yukawa couplings of the SO (12) model will be computed. 

To proceed let us consider the zero mode equations corrected up to 0(e 2 ) and in the 
absence of worldvolume fluxes. Using (3.8) one finds that the F-term (4.21) reduces to 

dfhVpxy + im 2 q^(p)a P rn + 2e9 2 m 4 [t^dya^ - t y (p)d x a pfh ] = (5.38) 

where we have defined the quantities 

t m (p) = xs mx (p) + ys my (p) ; m = x,y (5.39) 

These t m (p) can be regarded as the components of a 1-form. The values of s mn = s nm 
given in table 1 imply that dt(p) = 0. From (4.22) we simply obtain d x a py = d y a px . The 
D-term equation is given by 

d x a px + d y a py - im 2 q^(p)(p pxy = (5.40) 

In the following we present the solutions in the relevant sectors, again dropping subscripts 
xy in ip pxy . 

Sector a + 

We obtain a localized solution with a a + y = and 

a a+x = --^e x ^ 2 H a (y,x,x) ; <p a+ = ( H a + ^d x H a ) (5.41) 

m z \ wrx J 

The constant A a is such that the D-term equation (5.40) is satisfied for any H a . We 
further demand that tp a + — > m*fi(y) when x — > 0. Substituting into the F-term (5.38) 
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and solving to first order in e we obtain 

i\ 2 a m 2 
3(A^ + m 



H a = f t (y) + 2e6 2 [ xa - -^— f.(y) ) (5.42) 



The F-term equation is satisfied to (9(e) for any A a . Hence, we expect the A a dependence 
to drop out completely in the computation of holomorphic Yukawa couplings. 



Sector b + 

There are localized zero modes with a b + s = and 
iX b m 



a b+y — o~ 



e x ^H b (x, y, y) ; <p b+ = ( H b + ^d y H b ) (5.43) 

V m y J 

Imposing that the D-term equation (5.40) is satisfied for any H b determines X b . It is also 
required that (p b + — > m*gj(x) when y — > 0. Inserting into the F-term (5.38) and solving 
to first order in e leads to 

yfi- m i + m^ {x) ) (5 ' 44) 

iX b m 6 , 

V9i\x) 



12(A2 + m 4 )^^ 
The F-term equation is verified for generic X b . 

Sector c+ 

The wavefunctions in this sector have the form 

iX r m 



e A c |x-»| Hc {x,x,y,y) (5.45) 



\r\x-v\ 2 ( TT , ^ C ^xH c — dyH c \ 

¥c+ = m*e Ac|x y| [H c + — = ^ — 

\ m 4 x — y j 

Requiring that the D-term (5.40) is satisfied for any H c fixes A c . When x — > y it must be 
that (p c + — > m*7 c . Solving the F-term (5.38) to first order in e we find 

H c = 7c [1 + 2e9 2 (x - y)u] (5.46) 
iX^m 2 iX r m 2 



v = 



8(A2 + m 4 ) 16(A2 + m 4 
The F-term equation leaves A c free. 
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6 Yukawa couplings 



Given the non-perturbative zero mode wavef unctions one may proceed to compute the 
non-perturbative corrections to the Yukawa couplings. For this one must consider the 
cubic terms in fluctuations that arise from the full superpotential (4.1), insert the values 
for the background fields and zero modes and, finally, compute the overlapping integrals 
that contribute up to 0{e). 

As before, we consider the case where, from all the possible holomorphic functions 9 n 
present in (4.1), only 9q and 82 are non-vanishing. It is then clear that the cubic terms 
arising from this superpotential have the structure 

y^cubic = yj/xubic + yj/cubic 

where M/™ 1310 arises from the tree- level term (2.1) and W 7 ^ 111310 from the sum of the two 
terms in (4.1) proportional to 9q and 6 2 . 

Let us analyze these cubic terms in some detail. We have that the tree-level coupling 
reads 

jycubic = _ im 4 I Tr (</? A a A a) (6.2) 
Js 

-Ho) 

and leads to the expression (3.32). If we just insert the tree-level wavefunctions ip 

into (3.32) we recover the tree-level Yukawa couplings that we have already computed in 

section 3.4, whose single non-zero coupling is given by (3.36). If we instead insert the 

corrected wavefunctions (5.1) there will be in addition two sets of 0(e) integrals, that 

-Hi) 

will contain two tree-level zero modes and one 0(e) correction 6n , with n = 0,2. The 
Yukawa structure that arises from (6.2) is then 

^0 = Y< 0) + Y$ + Y$ + 0(e 2 ) (6.3) 

where is the tree-level Yukawa and Yq^ an 0(e) correction linear in ^ p \ n - 

The corrected couplings induced by W\ depend explicitly on the background and read 



e / 9 2 STr {{$ xy )tp xy Da A Da + ip 2 xy Da A (F)) (6.4) 
Js 



where D m a = d m a — i[(A m ),a\. The non-perturbative piece depending on 8q does not 
yield a term trilinear in fluctuations, as explained in detail in appendix D. The Yukawa 

43 



structure that arises from this term is simply 

Yi = Y} 0) + 0(e 2 ) (6.5) 

where (0) indicates that we must only insert tree- level wavefunctions into (6.4). 

In order to compute the couplings arising from W^ 1310 it is useful to derive a more 
explicit expression for the integral (6.4), by applying some specific features of the 5*0(12) 
model. Recall that quarks and leptons arise from the sectors a+ and b+ and come in 
families indexed by % and j respectively. They couple to a Higgs in the sector c + when we 
insert the zero modes ip p E p , p = a+ , &+, c + in the symmetrized traces of the superpotential 
(4.1). To simplify the notation in this section we drop the subscripts that identify the 
subsectors with different hypercharge. 

Let us now consider an specific 7-brane background. Namely, we will consider (& X y) 
given in (3.3) and a flux of the form (3.9), which we will write as (F) = $Qf- It is then 
useful to define the quantities 

O a +b+c+ = 6> 2 STr (($ xy )E a +E b +E c +) = d 2 m 2 (xd xa + b+c + + yd ya+b+c +) (6.6) 
K a +b+c+ = 2 STr (Q F E a +E b +E c +) = -9 2 (d xa + b + c + + d ya + b + c +) (6.7) 

where the d ma + b + c + are the traces 

d ma + b+c+ = STr (Q m E a +E b +E c +) ; m = x,y (6.8) 

with the generators in the fundamental representation. The correction due to W 7 ^" 1310 is 
then found to be 

( y i)«+6+ c + = 2m * 6 / °a+b+c+ Da\ + A Da> b+ ip c+xy + K a+b+c + <p l a+ X y<pl+ xy Da c + A $ 
^ s 

+ cyclic permutations in a + , b + , c + (6.9) 

It is straightforward to generalize this formula for a more general flux of the form (3.21). 
This will however not be necessary for our purposes, since in section 6.2 we will compute 
the corrected couplings for the case (F) = which already captures the correction to the 
holomorphic Yukawa couplings. In appendix D we will explain how the formula is applied 
when fluxes are present. 
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To summarize, we have that the corrected Yukawa couplings can be expressed as 



+ 0{e 2 ) (6.10) 



y* -(y^Y + (y w Y + \(y w Y +(Y W Y 

w c + - ) a+b+c+ + ^0,0 ) a+b+c+ + |y 0)2 ) a+b+c+ + i/i ) a+b+c+ 

where both Yq ^ and Yq 1 ^ come from (3.32) after inserting the corrected wavefunctions 
(5.1), and respectively extracting the zeroth and first order contributions in e. Similarly, 
Y^ is calculated substituting the uncorrected wavefunctions in (6.9). Notice that Y$ 
represent the 0(e) corrections to the Yukawas due to 9 , that is the corrections that we 
would obtain if 9 2 vanished, and which we compute in section 6.1. Similarly, the sum 
in brackets corresponds to the corrections due to 9 2 , which are computed in section 6.2. 
Finally we will consider that #0 is lineal in x and y and that 9 2 is constant, as this is enough 
to generate all the possible corrections to the Yukawa couplings below 0(e 2 ). With these 
assumptions the coupling Y 3 + b+c+ remains uncorrected to 0(e 2 ) and it is given by (3.36). 

6.1 Couplings due to 6q 

The corrected Yukawas implied by 9 are easy to compute by substituting the zero modes 
of section 5.1 in (3.32) and evaluating the integrals. The only non- vanishing couplings 
turn out to be 

2 6 

fTTTtl 

Y a 2 ? b+C+ = —T^ f a+b+c+ lhhc{9 0x + 9 0y ) (6.11a) 

2 fi 
flT 771 

y a % +C+ = — fa+ b +c + lhhc9oy (6.11b) 

Y£ b+C+ = ^ f a+b+c+ lhhc9 Qx (6.11c) 

To obtain these couplings we have used the functions Ai and Bj given in (5.13) and (5.18) 
respectively. In particular, the constants a and b are such that the couplings Y 13 and 
Y 31 turn out to be flux independent, up to normalization factors. For instance, before 
substituting the value of a we find that 

2 4 
67T Tfl 

Yah +C+ = fa+b+c+lhhc [(1 + 2Qm% x - £m% y + a ] (6.12) 

We see that taking a = rnl(( 2 9 0y — 2( a 9 0x ) indeed leads to (6.11c). In appendix D we 
show how the couplings (6.11) can also be deduced from a residue formula. 

To conclude, we have derived a set of Yukawa couplings that are flux independent up 
to the O(e ) normalization factors 7* similarly to the structure obtained in [22]. 
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6.2 Couplings due to 62 

In this case we need to substitute the wavefunctions of section 5.3 in (3.32), and also in 
(6.9). Since we have computed the zero modes with fluxes switched off, the corrections 
to the Yukawa couplings are simply given by 11 



l l+x 



( y i)a+&+c+ = ~ 2m * e / dvol s a +b+c+ {d x a l a+x d y a> h+ - - d y a l a+s d x a 3 h+ - - d x a l a+y - d y a\ 
J s 

+ d y a l a+ yd x al + -) tp c + xy + cyclic permutations in a + , b + , c + (6.13) 

The final outcome for the couplings follows by adding the different contributions that 
depend on 9 2 , as shown in (6.10). In the 5*0(12) model that we are discussing such sum 
simplifies by virtue of the properties 

% % 

d X a+b+c+ = -g/a+6+c+ 5 d ya + b + c + = —f a +b+c+ (6.14) 

Evaluating the various integrals involving the explicit wavefunctions gives the non-zero 
couplings 

Y* + \ +c+ = ^ f a+h+c+ lhhcO W (6.15a) 

Y^ b+C+ = - ^ fa +b+c+ lhhc02o (6.15b) 

where 620 = —i&2- Notice that these couplings are completely independent of the param- 
eters A a , A;,, and A c , which are not determined by the F-terms. 

We expect that with fluxes turned on the couplings will have the same structure (6.15), 
namely a flux independent core multiplied by flux-dependent normalization factors 7'. 
Had we computed the wavefunctions and Yukawa couplings for (F) ^ 0, we would have 
probably arrived to an expression similar to (6.12) in which the flux dependence in the 
core actually drops out after some parameters take their prescribed values. To sustain 
this assumption, in appendix D we compute the holomorphic couplings in presence of 
fluxes. To this purpose we will derive and apply a residue formula which is explicitly 
independent of fluxes. The couplings determined in this way completely agree with the 
results in (6.15). 



11 A parameter 9\ in the non-perturbative superpotential (4.1) gives rise to a similar coupling but with 
#a+6+ c + replaced by #iSTr (E a +E b +E c +) 1 which vanishes for SO (12) [22]. 
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7 Quark-lepton mass hierarchies and hypercharge flux 

Gathering the results of eqs. (3. 36), (6. 11), (6. 15) one obtains that the physical Yukawa 
mass matrix has a structure of the form 



Y D/L = 2ttV 2 7 c 



( 0(e 2 ) 0(e 2 ) eg-^ al f9 0x \ 

0(e 2 ) eg-^ b (9 0x + 6 0y ) 6^7^20 



(7.1) 



where the coefficients 



@0x — m2 @0x @0y — Tf^Q^y #20 — g m 2#20 (7-2) 

are adimensional, and we have defined the quotient of scales 

e= £)' = (2 ^ a (7 ' 3) 

The second expression for g can be derived in the type IIB orientifold limit of the 5*0(12) 
model. There, as explained in section 4, a — (m/m st ) 2 describes the intersection slope 
of the 7-branes in units of the string scale m~ t 2 = 2-ko.' . In this limit we also obtain the 
relation m A st = g s (2n) 3 m A , and combining both results the second equality of (7.3) follows. 

As already explained, the holomorphic Yukawa couplings are independent from fluxes, 
this dependence only appearing in the normalization factors 7* b c . In particular the de- 
pendence on hypercharge fluxes is the only possible source of distinction between D-quark 
and charged lepton Yukawas which are equal before this flux is turned on. In what follows 
we will be using the uncorrected normalization factors as computed in appendix A, since 
these corrections would only induce terms of order e 2 in the physical Yukawa couplings. 
The same happens with the mixing discussed in section 5.2, which only affects the physical 
Yukawa couplings at higher order. 

Our expressions for Yukawa couplings apply at the unification-string scale, presumably 
of order 10 16 GeV, so that in order to compare with experimental fermion masses one 
needs to run the data up to the unification scale. An updated two-loop analysis for this 
running within the MSSM has been performed in ref. [44] from which we will take the 
data below. Here we will only discuss fermion masses for charged leptons and D-quarks, 
which are the ones relevant for the SO (12) case studied here. Table 5 shows the relevant 
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fermion mass ratios evaluated at the unification scale for various values of tan/3 (the ratio 
of the two Higgs vevs in the MSSM). We also show for reference the Yukawa couplings 
of the r lepton and b and t quarks. Recall that m rb = Y Tjb Vcos(3, m t = Y t Vsm/3, 



with V = a/K 2 + Vd — GeV. As emphasized e.g. in ref. [44], the Yukawa couplings 



tan/3 


10 


38 


50 


m d /m s 


5.1 ±0.7 x 10~ 2 


5.1 ±0.7 x 10~ 2 


5.1 ±0.7 x 10~ 2 


m s /m b 


1.9 ±0.2 x 10~ 2 


1.7 ±0.2 x 10~ 2 


1.6 ±0.2 x 10~ 2 


me/m^ 


4.8 ± 0.2 x 10~ 3 


4.8 ±0.2 x 10~ 3 


4.8 ±0.2 x 10" 3 


m^/m T 


5.9 ±0.2 x 10~ 2 


5.4 ±0.2 x 10~ 2 


5.0 ±0.2 x 10~ 2 


m b /m T 


0.73 ±0.03 


0.73 ±0.03 


0.73 ± 0.04 


Y T 


0.070 ±0.003 


0.32 ±0.02 


0.51 ±0.04 


Y b 


0.051 ±0.002 


0.23 ±0.01 


0.37 ±0.02 


Y t 


0.48 ± 0.02 


0.49 ± 0.02 


0.51 ±0.04 



Table 5: Running mass ratios of leptons and D-quarks at the unification scale from ref. [44]. 
The Yukawa couplings Y T f>,t at the unification scale are also shown. 

obtained at the GUT scale seem to depart from the predictions from the minimal £77(5) 
GUT, which yield Y T ^ e = Y b:Sjd . This happens not only, as is well known, for the first two 
generations but also for the third for which one has a sizable departure from unification, 
see table 5. On the other hand, for large tan/3, which is going to be our case as we will see 
momentarily, there are additional large threshold corrections to Y b from the low-energy 
SUSY thresholds. In particular the leading such corrections in the MSSM are given by 



6Y b = 



91 ^ 3 tan/3 - ^^tan/3 



12tt 2 mi 



327T 2 m ( - 



(7.4) 



where M 3 ,m- b ,m^ are the gluino, sbottom and stop masses and //, A t are the Higgsino 
mass and top trilinear parameter. These corrections may easily be of order 20% (see 
e.g. [44], [45]) and the sign depends on the relative signs of the soft terms. So from table 
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5 and taking into account these corrections we will take for the third generation ratio 

Y 

= 1.37 ±0.1 ±0.2 (7.5) 

at the unification scale. The low energy threshold corrections may render b/r unification 
[44,45] but only for particular choices of parameters, particularly of signs. 

Since we do not know the corrections of order e 2 in the matrix (7.1) we will not 
attempt to describe the first generation masses but will only require that one of the 
eigenvalues should be much smaller than the other two. We will thus only try to describe 
the hierarchies between the third and second generation. From the table one gets for the 
mass ratios 

^ = 4.8- 6.1X10- 2 , — = (1.4-2.1) xl(T 2 (7.6) 
m T m b 

for tan/3 = 10 — 50. One then has for the (lepton/quark) ratio 

^#^3.3±1 (7.7) 

m s /m b 

We would like to see whether such hierarchies arise in our scheme. 
7.1 The third generation Yukawa couplings 

In order to compute the physical Yukawa couplings we need to compute the normalization 
factors Y a ,b,c- 111 principle this requires full knowledge of the matter wavefunctions along 
the matter curves. If however we assume that the relevant wavefunctions are localized 
close to the (unique) intersection point of the three matter curves, one can use the local 
convergent expressions for these normalization factors given in appendix A. As we said, 
we will be using the uncorrected normalization factors shown in that appendix, since the 
corrections would only induce terms of order e 2 in the physical Yukawa couplings. Then 
the leading contribution to the third generation Yukawa couplings is given by the 33 entry 
of the above matrix. Using (3.36) together with (A.65-A.67) we obtain 



. , ^ 2 (2A ai , 2 ±^ 2 (i + C 2 12 ))V /2 / ^ 3 (-2A fe2 , 3 ±gp 2 ' 3 (i±CL)) v " 

n,r = M ••'.'/•-) ' a I — d — M ■ ^ , ' I I — 

(2C C + q c P )(q c P + 2C - 2A C ) + (q c s + A C ) 2X ' 1 



^ 4 +a 2 1 , 2 (i±c\ 2 ) ) \ ^ 4 + ^, 3 ( l +c ft 2 2 , 3 : 



m 4 + C 2 + (Cc - Kf 



(7.8) 
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where the qp, A's and £'s are defined in appendix A and all fluxes, as in previous sections, 
have been taken constant in the vicinity of the intersection point. 

To estimate the value of the couplings we assume that the fluxes are such that g| ~ 0, 
q b s ~ 0, and q° s <C qp, as we will indeed find in our numerical fits. Using the results in 
appendix A, with M xy = 0, we then find the approximate result 

(ai, 2 b 2 , 3 c \ 1/2 
^^-) <"> 

The eigenvalues A are expected to be generically of the order of the charges qp ~ M. 

Taking q° s <C qp then shows that Y b / T is proportional to ag^J 2 with a small coefficient 

0(1). The intersection slope a is assumed to be small whereas g s is constrained as we 

now explain. Flux quantization requires {F) ~ 2tt, so that taking Vs 2 — V^ 2 implies 

M ~ Ny — Ny — (2tt)/V^ 2 . Next, the volume of the 7-brane surface S enters in the 

perturbative equality for the unification coupling (see e.g. [1]) 

27i 2 g s 
m%V s 

Setting a G ~ 1/24 leads to the estimate for the fluxes 



«G = TTTT ( 7 - 10 ) 



M (2a G \ 1/2 _ 0.29 



m st V 9s / g s ' 
Having diluted fluxes imposes M < m 2 st . We then conclude that g x J 2 cannot be arbitrarily 
small. The conditions of small intersection slope and fluxes are needed to justify the 
effective description of the 7-brane theory. 

Although (7.9) is just an approximation, it indicates that in the present scheme with 
the wavefunctions localized at the matter curve intersection point the third generation 
Yukawa couplings are large, of the same order of the Yukawa coupling of the top quark 
which we know on phenomenological grounds is of that order (see table 5). So e.g. in a 
MSSM scheme one expects a large tan/3 ~ m t /m b ~ 20 — 50. In the computations below 
we show that the above qualitative statements remain true in our F-theory scheme. 

7.2 Hierarchies of fermion masses. 

We have discussed above the Yukawa couplings for the third generation. The corrections 
discussed in previous chapters are in principle able to generate Yukawa couplings and 
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masses also for the first two generations. We would like to explore now to what extent the 
above results could be able to describe the observed structure of hierarchical lepton and 
D-quark masses. We will be interested now on the relative hierarchies among different 
D-quarks or different charged leptons, so that we will study the matrix 



Y 

y33 



la 



(7.12) 



where we have divided the matrix (7.1) by the largest F 33 entry. It is easy to check that 
this matrix has eigenvalues 

Ax = l + 0(e 2 ) 

A 2 = eg- ll ^\{9 Gx +~9 Qy ) + 0{e 2 ) 
A 3 = 0(e 2 ). 

This is interesting since we automatically get a hierarchy of masses of order (1, e, e 2 ) from 
the start, without any further assumption. Note that if we had 9 = and we were left 
only with the corrections from 9 2 , the matrix would be still rank one up to order e 2 and 
the non-perturbative corrections would be unable to create the desired hierarchies. So in 
order to obtain hierarchies it turns out to be crucial the presence of a non-constant 9 as 
studied in the previous sections. 

Identifying the first and second eigenvalues with the third and second generations one 
gets to leading order in e at the unification scale. 





( 


m b 






( q a P 2 qp 3 







\ 1/2 

j eg~\9 0x + 9 0y ) (7.13) 

\ 1/2 

J eg-\9 0x + 9 0y ) (7.14) 



where we have used 7 2 /7„ = {qp/ m V) 1 ^ 2 an d llhl = {~qp/ m T) 1 ^ • Using the expressions 
above together with the charges in table 2 we get for the ratio of ratios 

1/2 

(7.15) 



mJm T _ ( (M + \N Y )(M + N Y ) 



m s /m b y(M - \N Y ){M + ±N Y ) 
This ratio is interesting because the dependence on the non-perturbative correction and 
extra flux- dependent factors cancel out yielding a result which only depends on the ratio 

51 



m s /m b 
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Figure 2: at the unification scale as a function of Ny/M. The band shows the 

region consistent with the running of the low energy data up to the unification scale taken 
from ref. [44]. The dashed line is the Georgi-Jarlskog value [46]. 



Ny/M. Experimentally one has from table 5 that 
constraint 

~ ~ M 



m M /m T 
m a /m b 



3.3 ± 1 so that one has the 



1/2 



v 3M A 1 6M- 



3.3 ± 1 . 



(7.16) 



This condition is displayed in figure 2. One observes that agreement with experiment 
may be obtained with a ratio of fluxes Ny/M = 1.8 ± 0.6, independently from the value 
of the rest of the parameters. Note that conditions for consistent local chirality are 
M y + q Y + N Y > and M x + q Y + N Y < which are satisfied as long as -1 < N Y /M < 3. To 
reproduce the particular D-quark and charged lepton hierarchies we can substitute the 
last result in (7.13) and (7.14) which can be rewritten as 

1/2 



m s 




M 
mi 



M 
mi 



+ 0, 



Oy) 



7f).t: 



(7.17) 
(7.18) 
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From table 5 we see that — ^ = 5.4 ± 0.6 x 10 so that we get an estimate 

M 



rn 



2 e (6 0x + 6 0y ) ~ (2.3 ± 0.2) x 10" 2 (7.19) 



where we have used the central value N Y /M = 1.8 and the definition of g. Thus taking 
N Y /M ~ 1.8 and the non-perturbative correction of the size given by eq.(7.19) one obtains 
values consistent with the experimental m^/m T and m s /m b ratios. Note that the number 
in eq.(7.19) is quite small, consistent with a non-perturbative origin. 

7.3 b — r (non)-unification 

Hypercharge fluxes also violate the equality of the r and 5-quark Yukawas at unification. 
Indeed, using eq.(3.36) one gets 

' ,: " : (7.20) 



with the expressions for the 7's given in eqs.(A.65-A.67). We would like to see now whether 
one can obtain the result ^ = 1.37 ± 0.1 ± 0.2 discussed above for some choice of fluxes 
consistent with the equations for the Yukawa hierarchies of second to third generation 
discussed above. Recall that the latter require N Y /M ~ 1.2 — 2.4. We have as free flux 
parameters M, Ny, N a , since N b is determined by eq.(3.19) which sets Ny = 3(N a — Nj,). 
In figure 3 we show the value obtained for Y T /Y b as a function of the fluxes M and 
N Y and for values N a = ±1, which is sufficient to show the general behavior. This is 
done for N Y /M = 1.8 (upper plots) and 1.3 (lower plots), which correspond in turn to 
= 3, 2.2 respectively, yielding consistent 2nd to 3rd generation mass hierarchies. A 
first conclusion is that correctly yielding the latter hierarchy requires in turn Y T /Yb > 1, 
as observed. One can easily have flux choices with |f = 1.37 ± 0.1 ± 0.2, particularly for 
N a = -1 and N Y /M = 1.3 (lower left figure). 

An example of consistent parameter choices is (the fluxes are given in m 2 units) 

(M, N a , N Y , N Y ; e) = (2,0,0.1,3.6; 7.5 x 10" 4 ) (7.21) 

~ 1/2 ~ ~ 

where e = g s e (9 0x + 9 0y ). One obtains 

^, ^, Y b ) = (1.38, 1.7 x 10- 2 , 5.4 x 10" 2 , 0.66^V) (7.22) 
\Y b m b m T J 
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Figure 3: Y T /Yb as a function of M,Ny for N a = —1(1) left (right). Flux choices are consistent 
with = 3(2.2) for upper(lower) plots. 

in very good agreement with the experimental results in table 5. Note that the value 
for Yj, at the unification scale are consistent with those in table 5 for g s a 2 ~ 0.1 which 
suggest indeed a large value of tan/3 in the MSSM context. This is generally the case for 
all examples able to appropriately describe the mass ratios. Let us finally comment that 
the condition of diluted fluxes corresponds to e.g. M/m 2 st = (aM)/m 2 < 0(1), and the 
same for the rest of the fluxes. This may be achieved in the above example and also for 
the examples in figure 3 by considering an appropriately small value for a. 

One can repeat the analysis in this chapter for the case of the non-SUSY Standard 
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Model remaining below the string scale. Indeed, although we made our discussion in 
terms of superpotentials, the Yukawa coupling sector remains essentially unchanged in 
the presence of terms breaking SUSY below the string scale. In the case of the SM, 
extrapolating the low-energy masses up to the unification scale the ratios of the second 
to third generation masses remain similar to those shown in eq.(7.6), see e.g. ref. [47]. 
Concerning the b/r ratio one gets around the unification scale Y T /Y b = 1.73, with no 
relevant low-energy thresholds giving additional contributions. The flux analysis above 
would equally apply to this non-SUSY case and there are wide ranges of fluxes consistent 
with the fermion hierarchies and r/b ratio. The required fluxes tend to be however more 
diluted in this non-SUSY case. 

We conclude that instanton effects are able to generate the observed second to third 
generation hierarchies of charged leptons and D-quarks, via the superpotential deforma- 
tion (1.1). The relevant non-perturbative correction corresponds to the term proportional 
to 9 in (1.1), while the correction from 9 2 only affects the first generation. The hierarchies 
between second and third generations can then be easily understood in this scheme. One 
can reproduce the values (Y^/Y T )/(Y S /Y b ) = 3.3 ± 1.0 at the unification scale, consistent 
with the low-energy data and, at the same time, the b/r ratio with Yukawa couplings 
Y T /Y b (m st ) ~ 1.37±0.1 ±0.2, as obtained from the RGE in the MSSM. These two attrac- 
tive features are purely due to the hypercharge flux which explicitly breaks the underlying 
SU (5) symmetry, which otherwise predicts equal masses for D-quarks and leptons of each 
generation at the unification scale. The Yukawa couplings of the third generation are 
large, corresponding to large values of tan/3 (for not too small g s <J 2 ) within the context 
of the MSSM. In order to compute the masses for the first generation we would need to 
know the corrections to the Yukawa couplings at order e 2 which, although feasible, is a 
more involved task. 

8 Conclusions and outlook 

In this paper we have analyzed the generation and structure of Yukawa couplings in 
local F-theory models of 577(5) unification, taking into account the presence of non- 
perturbative effects. More precisely, we have shown how non-perturbative effects can 
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induce non-trivial Yukawas for the two lighter families of quarks and leptons, and naturally 
reproduce a hierarchical mass structure among the three families. In this sense, our work 
can be thought of as a continuation of [22], in which the same kind of approach was 
applied to a toy model of unification. While most of our previous conclusions remain 
true, there are several important new features that are crucial to understand the physics 
of Yukawa couplings in the present context of SU (5) unification. 

First, in a realistic F-theory GUT model the non-perturbative effects on 7-brane 
Yukawa couplings are non-trivially constrained. Recall that in [21, 22] one could un- 
derstand the effect of distant instantons on Yukawa couplings by simply adding the term 
W np = e J #iTr($F A F) to the 7-brane tree-level superpotential. However, such term 
identically vanishes when evaluated at Yukawa points of SO(2N) or E 6 j j8 enhancement, 
which are precisely the cases of interest for F-theory GUT models with chiral matter lo- 
calized at matter curves. As pointed out in [22] this does not mean that non-perturbative 
effects are trivial, but rather that one needs to consider the more general superpotential 
(1.1) and repeat the analysis with those terms that do not vanish. In the scheme of [21,22] 
the less suppressed of such terms is e J # 2 Tr($ 2 F A F), which we have indeed included in 
our analysis. We have however seen that the hierarchical structure that this term gener- 
ates for the masses of the three fermion families is of the form (0(e 2 ), 0(e 2 ), 0(1)), which 
is rather unsatisfactory from a phenomenological viewpoint. 

A more interesting result is obtained when the term e J 9 Tr(F A F) is included in the 
analysis. This term is in fact the least suppressed of the series (1.1) but becomes irrelevant 
for $o constant, which was the case for the scenario analyzed in [21,22]. However, in the 
context of realistic F-theory GUT models we may have a non-constant holomorphic 9q. 
Indeed, we have analyzed an SU (5) F-theory model in the vicinity of a Yukawa point of 
SO (12) enhancement, where the couplings 10 x 5 x 5 are generated. This model admits 
a natural type IIB realization, and there one can argue that 0(1) D3-brane instantons 
can induce a non-constant 6 . We have analyzed the 10 x 5 x 5 Yukawa structure with 
W np = e/ #oTr(F A F) + e / # 2 Tr($ 2 F A F), and obtained a fermion mass hierarchy of 
the form (0(e 2 ), 0(e), 0(1)) in a much better agreement with observation. It would be 
interesting to extend this analysis to Yukawa points of E Q j t $ enhancement, for which we 
expect a similar result. 
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Second, in the model of 5*0(12) enhancement one can analyze in detail the effect of 
the hypercharge flux Fy in the Yukawa couplings. As emphasized throughout the text 
the hypercharge flux is the only ingredient that breaks the SU(5) gauge symmetry of 
the model and, as a result, it should be able to explain the observed differences between, 
e.g., D-quarks and charged lepton masses. One can indeed check explicitly that due to 
the presence of F Y , all the SM particles have different internal wavef unctions. Hence, in 
principle the same SU(5) coupling 10 x 5 x 5 can lead to different SM Yukawas. It is 
however not clear that such differences in the wavefunctions is sufficient to understand 
the experimental data, as Yukawa couplings turn out to depend on the hypercharge flux 
in a rather subtle manner. 

Indeed, it was shown in [18] that holomorphic Yukawa couplings do not depend on 
worldvolume fluxes, including the hypercharge flux. In [22] this result was verified in the 
presence of non-perturbative corrections of the sort W np — e J 6 , iTr($F A F), with the 
help of the non-commutative formalism of [18]. In this work we have extended the above 
results and shown that the same statement is true for the non-perturbative correction 
W np = e j 6 Ty(F AF) + ej # 2 Tr($ 2 F A F) relevant for F-theory GUTs. Unlike before, 
this time there is no obvious interpretation of + W np in terms of non- commutative 
geometry, but we have nevertheless derived a residue formula for the holomorphic Yukawas 
in which the worldvolume flux independence is manifest. It would be interesting to study 
whether there is also a non-commutative structure behind W tTee + W np in this case. 

As a result, at the holomorphic level and within each family, all the Yukawas that 
arise from 10 x 5 x 5 have the same value. The only way in which one can reproduce the 
pattern of D-quarks and charged lepton masses observed in nature is via the wavefunction 
normalization factors that take us from holomorphic to physical Yukawa couplings. We 
have computed such normalisation factors for the local 50(12) model under study up 
to 0(e 2 ) corrections. Such computation is reliable when the matter wavefunctions are 
sufficiently peaked around the Yukawa point, which is a typical requirement for having 
reasonable Yukawas for the heaviest family. In fact, our analysis suggests that the Yukawa 
of the bottom is of the same order of magnitude as the top-quark Yukawa, which points to 
a large tan/3 within the context of the MSSM. Regarding the mass hierarchies within the 
third and second family, thanks to the hypercharge flux dependence on the normalisation 
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factors one is able to reproduce reasonable ratios (7.5) and (7.7) at the unification scale. 
This is a clear improvement compared to the classical 4d field theory models of 577(5) 
unification, in which such hierarchies can only be achieved by means of a complicated 
Higgs system and substantial threshold corrections. 
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A SO(12) wavefunctions at tree-level 

In this appendix we solve in for the zero mode wavefunctions of the SO (12) model of the 
main text. In particular, we will solve for the zero modes of the tree-level superpotential 
(2.1) and D-term (2.2). First we will find such zero modes in the holomorphic gauge 
(3.25), and show that they can be expressed as (3.30). We will then express them in a 
physical real gauge and find the normalization factors 7* p , Tjj and 7 Cr that correspond to 
canonical kinetic terms in the 4d effective action. 



A.l Zero modes in the holomorphic gauge 

As discussed in [22] and in the main text, in the absence of non-perturbative effects the 
7-brane bosonic (a^, ( -p xy ) and fermionic (ipm, Xxy) zero modes satisfy a Dirac-like equation 

D A ^ = (A.l) 

where \I/ is a 4-vector of wavefunctions and D A a 4 x 4 matrix of operators, both defined 
in (3.28). Massive modes satisfy instead 



DA f D A f 



(A.2) 



where Da^ is given by 



D A t 



V 






D- x 


Dy 


D- z 


D- x 





-D z 


Dy 


Dy 


D z 





-D 


D- z 




D x 






(A.3) 



/ 



In practice, in order to solve for the zero modes it is useful to consider eq.(A.2) with 



m p = 0. That is because the operator Da*Da takes the form 

D A f D A = -AI4 + M 



(A.4) 



where A is a Laplacian on the GUT 4-cycle S, and M is a matrix whose coefficients 
are 7-brane worldvolume fluxes and intersection slopes. For constant fluxes M and A 
commute, and so the first step to solve for the eigenmodes (A.2) is to diagonalize M. In 
the following we will solve for the 50(12) model zero modes for each sector of table 2. 
We will be sketchy, and we refer to [22,48] for more details on this class of computations. 
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Sector a 

For the sectors the matrix M in (A. 4) takes the form 



M„ 



± 



f M xy ^ 



— q a P q a s —im 2 

q a s q a P 
im 2 



V 



± 



M, 



xy 

m. 



where the entries q%, q a P take the values 



M xy I 3 







q a P 




-N a - iNy 


M-\N Y 


af 


-N a + \N Y 


M+\N Y 



(A.5) 



(A.6) 



and we have defined M and M xy as in (3.24). This flux matrix is diagonalized by the 
following unitary matrix 

I 1 o \ 

\\A X -\%Ai -i\%A* 
-A?CMi KCZM i\%QAz 
y im 2 Ai —im 2 A2 m 2 A 3 J 

where A" are the three eigenvalues of m a +, chosen so that A" < A^ < A3, and satisfying 
the cubic equation 12 

(A^) 3 - 2M xy (Xf) 2 - [m 4 + (q P ) 2 + (q*) 2 - M 2 y }\? + m\q P + M xy ) = (A.7) 
Finally, we have defined A { = ((A?) 2 + (A^Q 1 ) 2 + m 4 )~ 1/2 and 



C = 



(A.8) 



\?-q a p - M xy 

It is easy to see that the lowest eigenvalue of M for the sectors a+ is given by A" — M xy , 
and by M xy — A3 for the sectors a~ . The potential fermionic zero modes for any of these 



12 In terms of the quantities A + , A_ defined below table 3 we have that A + (a+) = A" and A_(a p ) = A3. 
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sectors should then be proportional to the corresponding eigenvectors of M. Indeed, 
considering the sectors a+ , the Ansatz for a fermionic zero mode would be 

/ \ ( n \ 



1 






m 2 



(\l-q a p -M xy )m? 



(A.9) 



v°y v 1 / 

where y' + and y + are scalar wavefunctions that differ by a constant factor. With this 
Ansatz eq.(A.I) translates to 







(A.IO) 



where the differential operators D a correspond to the entries of the Dirac operator T>x = 
Ja p D A J ap in the basis rotated by 3 ap . Namely, we have that 

( \ ( A-y [\\D X - XiCiDy + im 2 D z \ ^ 

A 2 [-\ a 2 D x + X-QD y - im 2 D z ] (A.ll) 
v A 3 \-i\%D x + i\lQD y + m 2 D z ] j 



V 



V J 



f A 1 [\ a 1 D x -\ a l (ZD y -im 2 D z ] N 

„4 2 [-AfAg + A^A/ + im 2 D 2 ] 
v A 3 [iX a 3 D x -i\ a 3 QDy + m 2 D,} } 



(A.12) 



The general solution to eqs.(A.lO) in the holomorphic gauge (3.25) is then given by 



a P \ \ a — n a — 



-X 



(A.13) 



K-q a P - M xy 

where f a +, p = 1, 2 are holomorphic functions to be determined by boundary conditions. 
Taking \ ap = A" and Ca p = Ci we obtain the zero mode wavefunctions 

/ n \ 



hoi 







a p m 2 



\ 



e^-^f at (y + ( ap x)E ai 



(A.14) 



/ 



that is indeed of the form (3.30), with A + = \ ap . 
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Sector b 



For the sectors bf, q = 1,2, 3, the flux matrix M takes the form 



( —M xy ^ 



M, 



■q b P Is 



\ 



q b P im 2 






-im 2 
where q b s , q b P take the following values 



—M. 



xy 







m b + + M xy I 3 





4 




bf 


N b + §iV y 


-M + \Ny 


bf 


N b - lN Y 


-M - \N Y 


bf 


N b -N Y 


-M-N Y 



(A.15) 



It is diagonalized by the unitary matrix 

d 



\ 



(A.16) 





\\Ai -\ 2 A 2 -i\%A 3 

-AjCMi KQA 2 i\ a 3 QA 3 

y im 2 A\ -im 2 A2 m 2 A 3 J 

where \\ < \ b 2 < A| are the eigenvalues of m b + and so they satisfy the equation 

(A*) 3 + 2M xy (X b ) 2 - [m 4 + (q b P ) 2 + {q b s ) 2 - M 2 xy }X b - m\q b P + M xy ) = (A.17) 
Finally, we have defined A l = ((A-) 2 + (A?Cf) 2 + m 4 )~ 1/2 and 



c b = 



Qs 



\ b + q b P + M xy 

As before, the appropriate Ansatz for fermionic zero modes in the sectors bf is 



(A.18) 



%t = 3b 



1 



V / 
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Xb+ E b+ 



(A.19) 



Therefore, the wavefunctions Xb+ mus t satisfy the differential equations 

V xX b + = VyX b + = V- zXb + = 
where V A = D A J& 9 is the Dirac operator in the rotated basis: 

( V x \ I A, [X\D X - \ b ( b D y + im 2 D z ] ^ 



(A.20) 



Vy 



(A.21) 



(A.22) 



A 2 [-\ b 2 D x + \ b 2 ( b D y -im 2 D z ] 
\ A 3 [-i\\D x + i\ b ( b D y + m 2 D z ] J 

' A l [\\D x -\ b l C b l Dy-im 2 D z \ \ 

A 2 [-\ b 2 D x + \ b ( b 2 D y + im 2 D z ] 
K A 3 [i\ b D x - i\ b ( b D y + m 2 D z ] J 
It is easy to check that the zero mode equations in the sector 6+ are identical to the sector 
a+ provided that the following replacements are performed 

x y A" <->■ \\ Q Cf 
9s ^ 9s 1p ~<l b p M xy ~M xy 
Thus, the solution of (A.20) in the holomorphic gauge is 



(A.23) 



4 



9s 



-y 



(A.24) 



\ b + q b P + M xy 

where f b +, q = 1,2,3 are arbitrary holomorphic functions. Taking A^ = Aj and = ^ 
we obtain the zero mode wavefunctions 



vl/ 11 ? 1 = 





-Ck ^ 

iX bq 



\ 



\ 



1 



e \y(y-(^) f {x + (bq y )E 



(A.25) 



/ 



that is indeed of the form (3.30), with A + = \b q - 



Sector c 



For this sector the flux matrix reads 

/ 



M c ± = ± 









-Qp 



9s 
9p 



-im 2 im 2 



im 
—im 2 




\ 



± 




m - 



(A.26) 
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with q s , q c p now taking the following values 





fs 


Qp 


of 


N a -N b - ±Ny 


-iJVy 


4 


N a -N b + \N Y 


\Ny 



The unitary matrix diagonalizing M c ± is 

(l 



V 





o CMi Cl^2 Cl^2 

(Cr + AJ)A (Q + \ C 2 )A2 (Q + K)^ 

im 2 A\ im 2 ^ 2 im 2 Az J 
with AJ < X c 2 < A3 being the solutions to 

( A c)3 _ [ 2m 4 + ((? c } 2 + (? c } 2] A c + 2m 4 g c = 

and where we have defined A c a 



(A.27) 



((Cf) 2 + (A<-C) 2 + ™T^and 



A?(A< 



-1/2 



m - q c s) 



(A.28) 



(A.29) 



A? + qp-q c s 
In the last equality we used (A.28). 

As discussed in appendix B, for q c s 7^ this sector yields a chiral spectrum. Let us 

focus on the case where q c s > and look for zero mode solutions in the sector c+. The 

appropriate Ansatz is 

/o\ 



1 


V / 



X c +E c 



(A.30) 



with the wavefunctions x c + satisfying 

V xXc+ = V vXc+ = V zXc+ = 

and where V& is the Dirac operator in the new basis, namely 

( V r .\ I Ax {-CtD x + (XI - Cl)D y + im 2 D z ) ^ 
A 2 (-QD X + (X c 2 - Q)D y + im 2 D z ) 
K A 3 (-QD X + (XI - Q)D y + im 2 D z ) } 
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(A.31) 



(A.32) 



Vy 



' A 1 {-CID- X + X1- (Ci) Dy - tm 2 D,) \ 

A 2 + (A c 2 - Q)D y - tm 2 D- z ) 

\ A 3 {-QD- X + (A3 - Q)D y - tm 2 D 2 ) J 



(A.33) 



The solution to (A. 31) then reads 

= e l*-v)(&-M-CM) f ((A c _ C c )x + C c y) 



(A.34) 



with f c + holomorphic. Again, this solution fits into the general expression (3.30). Taking 
XI = \ Cq and ( Cq = (1 we find the solution 

/ n \ 



hoi 





m 2 



e (*-y){^-^ q -t Cq )y) f4 ( (Ac? _ ^ + (cq ^ E (A _ 35) 



X,Y bosons 

Let us finally consider the sectors X ± , Y ± in table 2, arising from the breaking SU(5) — > 
SU(3) x SU(2) x U(1)y triggered by the presence of the hypercharge flux. Such hyper- 
charge flux is typically chosen so that no zero modes arise from this sector, as they would 
correspond to SU (5) exotics. However, there will be massive modes, and in particular a 
tower of 4d massive gauge bosons that can be identified with the X ± , Y ± bosons of 4d 
517(5) GUTs. 

Unlike the cases analyzed previously, the sectors X ± ) Y ± do not correspond to any 
matter curve, and so they contain wavefunctions of the bulk of the 4-cycle S. As a result 
the quantity m 2 does not appear anywhere in the flux matrix, which reads 



/ 



M 



X,Y± 



= ±- 



V 





-N Y N y 

N Y N Y 





\ 



(A.36) 



/ 



and that as expected becomes trivial for vanishing hypercharge fluxes Ny and Ny- This 
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matrix is diagonalized by 



where 



iVv 



c = 



s = 



10 

c -s 

s c 

1 

N Y - A 



\ 



(A.37) 



A = y Ny + N Y (A.38) 



N Y + (N Y - A) 2 ^N Y + (N Y - A) 2 

It is easy to obtain the lowest eigenfunction of the scalar Laplacian A. Indeed, computing 
the rotated Dirac operator V A = y D A Jx,y gives us 



( V„ \ ( cD x + sD y ^ 



Vy 

V 



\ 



-sD x + cD y 



(A.39) 



/ 



/ T>- \ ( cD x + SDy \ 



V,-, 



\V- Z J 



\ 



(A.40) 



/ 



—sD x + CDy 

D- z 

Following appendix A of [22] one can convince oneself that the eigenfunctions of the 
Laplacian with lowest eigenvalue satisfy 



V x Xx,y+ = V yXx,y+ = V s Xx,y+ = 

VxXx,Y- = VyXx,Y- = VzXx,Y- = 

Notice that if we define the rotated holomorphic coordinates u and v as 



(A.41) 



u = cx + sy 



v = —sx + cy 



(A.42) 



then in the holomorphic gauge the above operators read 

5 
6 

5 



V x = d u ± | Att V x = d u 



(A.43) 



Vy = 8 V =F | At) Vy = 

for the sectors X, Y^ 1 respectively. Hence, the solution to (A.41) in the holomorphic gauge 
is given by 



AM 2 



hoi 



Xx,Y" = e 6 /X,Y- l«, Uj 



(A.44) 
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with f x Y ± depending on one holomorphic and one anti-holomorphic variable. Finally, it 
is easy to see that the lightest 4d gauge boson of this sector 

/ 1 \ 






V / 



has (A. 44) as scalar wavefunction and satisfies the equation (A. 2) with 



(A.45) 



(A.46) 



As N Y , N Y are flux densities over the 4-cycle S, we typically have m X y ~ Vol(S') 1//4 . 



A. 2 Zero modes in the real gauge and normalization 

While the holomorphic gauge simplifies the computation of wavefunctions, it is not obvious 
to see from (A. 13) or (A. 14) if these zero modes actually exist. In this sense taking a 
real gauge provides more information, since one can study the local convergence of the 
wavefunction [22] and describe a set of wavefunctions that is locally chiral [36]. In the 
following we will compute the zero mode wavefunctions in the real gauge for the sectors 
discussed above, and compute the normalization factors for each of them. 
Instead of (3.25) let us consider the real gauge 



(A) 



real 



| [x(Q P - M xy Q F ) - yQs] dx - \ [y{Q P + M xy Q F ) + xQ s ] dy 
[x(Q P - M xy Q F ) - yQ s ] dx + l [y(Q P + M xy Q F ) + xQ s ] dy 



"2 



(A.47) 



= (A) hol + dQ 



where 



\x\ 



Q P + (xy + yx) Q s + M xy (\y\ 2 + \x\ 2 ) Q F \ 



(M x \x\ 2 + M y \y\ 2 ) Q F - N Y (\x\ 2 - \y\ 2 ) Q Y + (xy + yx) Q s 



(A.48) 



As discussed in [20], one can relate the wavefunctions in the holomorphic gauge to the 
ones in the real gauge by the simple formula 



^rcal ^iQ-qjhol 
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(A.49) 



where Qs, Qp and Qf act on the roots E p as in (3.26) and (3.11), respectively. This only 
changes the scalar wavefunctions, and so in the real gauge our zero mode solution for the 
sectors a p still reads (A. 9) but now with x a + & ven by 

^real = e _ M*-y* Y |x|a _ M y+q yN Y ]y? _ q% Re{xg)eXapX{ x_ Capg) ^ ^ + ^ ^ ^ 

Notice that the exponential prefactor of the wavefunction allows to compensate the di- 
vergent behavior of the holomorphic function f a + whenever \ ap — M x + qyNy < and 
My + qyNy > 0. One can check that the first condition is automatic, while the second 
needs to be imposed in order to obtain locally convergent zero modes in the sectors a+ . 
The same condition M y + qyNy > forbids the presence of zero modes in the conjugate 
sector a~ , at least at this local level. 13 Similarly, for the sectors 6+ we have 



and so for having locally convergent zero modes we need to impose M x + qyNy < 0. 
Altogether we have that M x + qyNy < < M y + qyNy, yields zero modes in the sectors 
dp and bg, as assumed in the main text. For the remaining sectors we have 

^real = g ( \N Y He) N 2 +(A C -C C - § N Y q Y ) \yf- ( _ Cc+ A c )xg- ( \q%Hc)xy j ^ + _Cc_ y ^ 

= e-^\ 2 ^)f^ Y+{ u,v) (A.52) 
X^- = e-A^(H 2 +M 2 ) /xiY _ ( U ,v) 

Following section 3.4, let us take the basis of holomorphic functions f p (v p ) = Y P m i~ lv p _i 
for the zero modes in p = a p ,6+ . Here i — 1,2,3 labels the three different families and 
7* is a normalization factor that we want to fix in order to have canonical kinetic terms 
for our matter fields. The normalization condition is given by 

— >real — >-real f — >-rcal — y t real 

i^ pi \^ pj ) = ™l Tr(^ pi • ^ p . ) dvol 5 (A.53) 

J s 

= 2ml0 p \\ 2 I X P {x j p y^o\ s = 5.J 
J s 



13 Our ultra-local description is insensitive to global features of the GUT 4-cycle S, and in principle 
there could be zero modes in the sectors that locally seem divergent but are globally convergent, and 
the other way round, see [36] for a discussion of this point. In the following we will assume that the three 
zero modes that represent the three families of 5's and 10's are local zero modes in the language of [36]. 
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where we have used that in our conventions Tr (EpE^) = 25 pa , and defined the vector 
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m 2 






i\y 






m 2 
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/ 



(A.54) 



whose entries are given in table 3 or in eq.(3.33). 

Notice that the normalization condition is automatically satisfied for i ^ j. Indeed, 
all the family dependence of our matter fields is encoded in the scalar wavefunction x% 
which in a real gauge can be written in the form 



x * p = 7 Jm*- i e-( a ' x l +%l +cx y +dyS \ax + (3y) 3 ~ l i = 1, 2, 3 



(A.55) 



for a certain set of constants a, b, c, d, a, (3 that vary for each sector. One then has that 
the norm of this scalar wavefunction 



= ml x;(4)*dvolc 
Js 



(A.56) 



vanishes identically for i ^ j, as the integrand is not invariant under the U(l) action 
(x,y) — > e ta {x^y). For % = j we can compute this norm by extending the integration 
domain from our local patch to C 2 

||V|| 2 = W\\mlf- 1 [ e -(2«W 2 +2%P + 2 (c+ d)Rc^]) (| ax + /3y | 2 )3-Mvol c2 (A.57) 
Jc 2 
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One then obtains 14 

NXap " " m^lr^lMS-O! \ "': J 1 ' 

IIXbJI = ^kims-o! 1 — ~J (A ' 59) 

llXcir 2 = -^^^((2C c + iV Y g Y )(iV y g y + 2C c -2A c ) + (g5 + A c ) 2 )(A.60) 

where, as above, for the sectors p = c r , X, we have taken the holomorphic function f p 
to be a constant. On the other hand we have that 

||^ a || 2 = m- 4 (m 4 + A 2 (l + C 2 )) (A.62) 
||^ b || 2 = m~ 4 (m 4 + A 2 (l + C b 2 )) (A.63) 
||^ c || 2 = m" 4 (m 4 + C + (C c - A) 2 ) (A.64) 

One finally finds the following normalization factors 

4 _o./n\ , _a/i , a2\N / _n, \ 3-i 



It: 



|2 



1 / m \ 4 (2C c + gfQ(gf> + 2C c -2A) + (g| + A) 2 
l7crl 2vr 2 VmJ m 4 + C 2 +(Cc-A) 2 1 ' 

where for simplicity we have set M xy = 0. 



14 We have made use of the following formulas arising from (A. 57) 
\K\\M\~ 2 = {i ^f,™l d?f ("V + b2 " 1 + «W C + d ? + ° i>c " 3 " 2 - ^ + d)aD S - 2i>(c + d)a'(l)) 



3 1 1 2 I 3 1 — 2 



11x^11^ 



(4a6 - (c + d)2) 
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B Fluxes in the SO (12) model and local chirality 

In section 3 we have made some choices in order to construct the SO (12) model analyzed 
in the rest of the paper. More precisely, given the set of matter curves that describe a 
point of 5*0(12) enhancement, we have chosen a worldvolume flux F can reproduce the 
matter spectrum of the MSSM. In this appendix we discuss and motivate the restrictions 
that have been imposed on such worldvolume flux, with some emphasis on the piece that 
corresponds to the hypercharge flux F Y which plays a particularly important role in the 
computation of physical Yukawa couplings. 

Recall that the chiral spectrum of an F-theory model depends on the integral of F over 
each matter curve E Q , while the choice (3.21) only describes F in a local patch around 
the point of 5*0(12) enhancement. Hence, in principle there could be no relation between 
(3.21) and the actual spectrum that each matter curve holds. However, there is a simple 
way to relate the two, namely by the concept of local chiral modes discussed in [36]. 
Indeed, following [36] let us define local zero modes as those whose internal wavefunction 
peak in a neighborhood of the point of 50(12) enhancement. We may then demand 
that the MSSM spectrum localized at the curves a, b and c arises from local zero modes. 
This latter assumption is in fact implicit in our discussion of zero modes, since otherwise 
computing the wavefunction normalization as in (A. 57) would not be accurate. 

The concept of local zero modes is important in our discussion because by assumption 
these modes are affected by the local density of worldvolume flux F at the point of 
50(12) enhancement. Hence, having the appropriate spectrum of chiral modes puts local 
restrictions on F. For instance, assuming local zero modes and following the computations 
of subsection A. 2, one arrives to the conclusion that the local flux densities in (3.9) must 
satisfy M x + qyN Y < < M y + qyN Y to have the proper spectrum on the curves E a and 
E b . In the following we would like to show how also having the correct local spectrum for 
the curve E c implies adding the extra pieces of worldvolume flux (3.12) and (3.13). 

Indeed, as pointed out in the text with only the flux (3.9) the curve c remains locally 
non-chiral, while in principle we would like to have a chiral spectrum that hosts the 
577(5) Higgs 5 and then, when including the hypercharge flux Fy, a chiral spectrum for 
the Higgs doublets and a non-chiral one for the triplets [5]. To have a chiral spectrum at 
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S c we simply need to add an extra piece to the previous worldvolume flux, and a natural 
candidate is the non-diagonal flux (F 2 ) given in (3.12) because adding this flux does not 
change the number of families in the curves a and b. This should be so because 



(F 2 )\ x =o = (F 2 )\ y=0 = 



(B.l) 



and so adding this flux does not change the integral of (F) over S a or E 6 . On the other 
hand, the pull-back of (F 2 ) on x = y does not vanish, and so turning on F 2 should have 
an effect on the curve c. 

A simple way to check how worldvolume fluxes affect the spectrum of each of the 
curves is by considering the operator 



xy 



yx 

F - 
r yy 

D $ 



-A 



xy 



-DyHy 




(B.2) 



/ 



y J-^x^xy ^y~*xy 

that arises when considering (A. 4), as seen in the previous appendix. Now, det m p is 
the product of the three eigenvalues of the matrix m p , which from [22] we know enter 
into the algebra of creation and annihilation operators for each curve. In particular, it 
follows from [22] that for a curve with a (locally) chiral spectrum we need the three of 
these eigenvalues to be non- vanishing, and so det m p ^ 0. 

In fact, one can check that det m p matches the definition of local chirality index of [36]. 
In the following we will analyze the value of this index for each sector of the S'0(12) model. 
From such analysis we will obtain a set of restrictions that must be imposed in the local 
worldvolume flux in order to have an acceptable local chiral spectrum. 



sector a + 



For this sector, in the presence of the fluxes (3.9) and (3.12) only we have that 



V 



M x -N a -im 2 


2 



-N a M y 



im 



\ 



J 



(B.3) 



Hence, before adding the hypercharge flux, we see that det m a + does not depend N a , 
and that det m a + = for M y = 0. The first observation means that, as expected, (F 2 ) 
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does not affect the chiral spectrum of this curve, and the second means that if we switch 
off the flux My then the curve a locally becomes a non-chiral sector. From subsection 
A. 2 we know that a positive number of chiral zero modes is obtained for M y > 0, which 
corresponds to det m a + < 0. 

When adding the hypercharge flux, this flux matrix becomes 



m a + 



/ 



V 



M x - q Y N Y 



Is 



im 



M y + q Y N Y 




—im 





\ 



/ 



with q Y (af) = — | and gy(a^~) = \. The local chirality index then reads 

detm a + = -m 4 (M y + q Y N Y ) 



(B.4) 



(B.5) 



and so it has the same sign as before as long as M y + q Y N Y > 0. Note that the flux 
density N Y that also enters into the hypercharge flux does not modify the local chiral 
index. However, from the computations in appendix A one can see that both flux densities 
N Y and N Y enter into the expression for the wavefunctions of this sector, the latter one 
through the quantity q$. 



sector b + 



With only (3.9) and (3.12) we have 



m b 4 
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— Mr.. 



N h 



N b -M y 
—im% 
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(B.6) 



Again, det m b + does not depend on the flux density N b , but rather on M x . Recall that 
obtaining local zero modes in this sector corresponds to taking M x < 0, which again 
translates into det m b + < 0. 

After adding the hypercharge flux we have 



m b 



( —M x — q Y N Y q b s 

-M y + q Y N Y im 2 
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—im 







(B.7) 
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where now qy(bt) = I, <2v(&2~) = ~~ | an d 9V(&3 ) = — 1- The local index reads 



detm b + = m 4 (M ;E + q Y N Y ) 



(B.8) 



and so we need M x + qyiVy < 0. This condition and the one for the sector a + are 
automatically satisfied whenever M x < < M y and \M X \, \M y \ > \N Y \. 



sector 



Finally, for this sector we have that 



m c + 



N n -N h im 2 



N a -N b 



—im 
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—im 



im 



(B.9) 
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when only the fluxes (3.9) and (3.12) are present. We see that 



detm c+ = -m 4 2(N a -N b ) 



(B.10) 



and so, for N a ^ N b we should have a chiral fermion localized in the curve c. As before, 
the condition to have a chiral spectrum of left-handed chiral multiplets arising from the 
sector c + (as opposed to arising from c - ) can be expressed as det m c + < 0, which by 
(B.10) is equivalent to iV a > N b . 

Once introduced the hypercharge flux (3.13), we have instead 
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-qyN Y q c s im 
qyNy —im 2 



\ 



—im 2 im 2 



(B.ll) 



and so the local chirality index becomes 



det m r 



-m*2q c s = -m 4 2(iV a -N b + q Y N Y ) 



(B-12) 



Since qy(cf) = — | and gy(c^) = |, one may achieve a locally non-chiral sector of Higgs 
triplets by simply imposing N Y = 3(N a — N b ). This constraint and the previous condition 
N a > N b imply that q c s > for the sector cj, which means that there is indeed a Higgs 
doublet in this sector. 
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C Non-perturbative superpotential and D-term 



The key ingredient that allows to obtain Yukawa matrices of rank higher than one is 
the deformed 7-brane superpotential (4.1) As pointed out in [21] and further discussed 
in [22,42,49], this superpotential arises once that one takes into account certain non- 
perturbative effects present in type IIB/F-theory compactifications. The purpose of this 
appendix is to recall the derivation of eq.(4.1), basically following [21,22], to which we 
refer for further details. In addition and based on the results of [42] we will show that the 
same class of non-perturbative effects leave unchanged the expression for the D-term. 

C.l Non-perturbative superpotential 

As reviewed in section 2, the basic setup of type IIB/F-theory GUT compactifications 
assumes that the GUT gauge degrees of freedom arise from a 4-cycle Sgutj a divisor of a 
compact threefold base B that is wrapped by a stack of 7-branes. This divisor intersects 
another set of divisors Si also wrapped by 7-branes, and it is at the intersection curves 
Sj where the chiral matter multiplets arise from. Near the Yukawa point where three 
of these matter curves meet one can understand this configuration in terms of a single 
stack of 7-branes wrapped on a divisor S and with a higher rank gauge gauge group like 
SO(12), E 6 , E 7 or E s . This last divisor and gauge group are the ones that enter into 
expression for the superpotential (2.1) and D-term (2.2), from which we can derive the 
zero mode equations of motion and the Yukawa couplings. The fact that a single 4-cycle S 
describes a set of intersecting divisors is encoded in the non-trivial profile of (<&). Finally, 
as pointed out in section 4, S and Scut need not be the same. In particular, in the 
type IIB description of the SO (12) model S should be associated with the divisor of an 
07-plane. 

In addition to all these divisors, there will be in general other set of divisors also 
wrapped by branes which may not be part of the MSSM sector, but instead the source of 
non-perturbative effects. Typical examples are 7-branes that give rise to a gauge hidden 
sector of the theory that undergoes a gaugino condensate, or Euclidean BPS 3-branes with 
the appropriate structure of zero modes to contribute the the superpotential of the 4d 
effective theory. In the following we will consider that one of such objects exists and that 



75 



is wrapping a divisor S np , and describe how the presence of S np affects the superpotential 
(2.1) for the gauge theory on S. 

Perhaps the most intuitive way to explain how the presence of S np modifies the 7-brane 
superpotential (2.1) is by considering the case of a gaugino condensate on S np . There we 
know that a non-perturbative superpotential of the form 

W np = fi 3 e-^ /r (C.l) 

is generated at the level of the 4d effective theory. Here / np is the gauge kinetic function 
of the stack of 7-branes on S np and r is the rank of their gauge group. Finally ji ~ 
is the UV scale at which / np is defined. This expression and all those below are also 
valid for the case where S np hosts an Euclidean 3-brane, by simply replacing / np by the 
complexified action of the instanton and setting r — 1. 

In general, / np will be a holomorphic function of the 4d chiral multiplets of the theory, 
which arise either from the bulk or from the 7-brane sectors of the compactification 

/np = T np + f£*°» (C.2) 

where the first term is the gauge kinetic function f 7np computed at tree-level, given by 
the complexified Kahler modulus T np = Vol (S np ) + i J s C 4 that corresponds to S np . The 
second term arises from threshold effects, and is a holomorphic function of the bulk/closed 
string fields and of the 7-brane fields. More precisely, if we consider the presence of a 7- 
brane in the divisor S we have that [21] 

/np" l0 ° P = -rlog-4 -^^STr(log/iFAF) (C.3) 

where A is a function of the bulk/closed string fields, and all the dependence of the 7- 
brane fields enters into the second term of this expression. There the key quantity is 
h, which is the holomorphic divisor function of the 4-cycle S np = {h = 0} sourcing the 
non-perturbative effect. While log/i is a scalar bulk quantity, when plugged into the 
expression (C.3) one should follow the prescription of [50] and consider its non-Abelian 
pull-back into S. That is 

log h = log h\ s + m- st 2 § xy [d z \og h] s + m; t 4 ® 2 x M\og h] s + . . . (C.4) 
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where as in the main text z is the holomorphic coordinate transverse to S, <& xy the 7-brane 
position field related to it and m~ t 2 = 2na'. Plugging this expression back to (C.3) we 
obtain 



1— loop 



■\ogA+ N m log ho + ^J2 ^ j} 9 ^ {h/h ° )]s STl (< ^ FAF) (C ' 5) 



where ho — f s h the mean value of h in S and Nd 3 = (87T 2 ) -1 J s Tr(F A F) e N the total 
D3-brane charge induced on S by the presence of the magnetic field F. Finally, inserting 
(C.5) into (C.l) we obtain 

-In 

I \f>Jn f r(h./h J ,'\]cTr(G. 

■ xy ■ 



W np = , 1 3 (^e- T » /r /i" m/ ')e X p 



m 



St 



87T 2 r J s 



[d z \og(h/ho)}sTr($ n xy FAF) 



(C.6) 



^e- T "< M/r f 1 + £ inr I (VMs STr F A F) + 



87r 2 r J 5 

Hence, up to a constant term we have that 

W np = mf-J s 6 n STr(^ y FAF) 

where we have defined 



(C.7) 



e = Ae' T ^ r h^ 3/r 

-nil 3 / 

g s fj, a /r 



n — 



[d> g (VMUo 



(C.8) 
(C.9) 



(27r) 2+ ^mt +2 " 

where we have used the relation m\ t = g s (27r) 3 mf derived in [22]. 15 Clearly, these expres- 
sions match the definitions of the text for the case r — 1. 



C.2 Non-perturbative D-term 

Interestingly, one can derive the same result (C.7) by replacing the presence of the non- 
perturbative sector in S np by a deformation of the compactification background, or in other 
words by 'back- reacting' the non-perturbative effect. This idea was put forward in [51] 
and further analyzed for gaugino condensate type IIB vacua in [42,49]. As discussed in 
[22] considering this deformed background simplifies the computation of non-perturbative 
7-brane Yukawa couplings. We will not follow this approach to rederive in detail the 
15 In fact, a factor of g s was missed in [22], which we incorporate here. 

77 



non-perturbative superpotential (4.1), referring the reader to [21,22] for the appropriate 
discussion, but will make use of it to derive the form of the D-term in the presence of the 
same non-perturbative effects. 

Indeed, in the scheme of [51] one may take the deformed background obtained from 
back-reacting a gaugino condensate and compute the D-term in such background. More 
precisely, in the language of generalized geometry, one needs to find the polyforms or pure 
spinors that define such background, and them derive the D-term (or F-term) equations 
from them. For the case of a gaugino condensate arising from a D7-brane, such pure 
spinor analysis was carried out in [42]. In particular, it was found that the pure spinor of 
interest to compute the D-term is given by 

e 2A lm = J+^AQ + ^(QAO) + ... (CIO) 
8 8 

which is basically eq.(5.15) of [42]. Here /3 is a real bivector 

(3 = /?( 2 >°) + /?(°> 2 ) = p^di A dj + c.c. (C.ll) 

whose (2,0)-component is holomorphic, and specifies the so-called /3-deformation of the 
background. In a given patch of the compactification (3 specifies (up to a constant) a 
holomorphic function Xo, since 

ip{l = ^m ijk dz k = ed Xo (C.12) 

Xo, in turn, specifies the deformations of the superpotential of the form (C.7). Indeed, we 
can define 

9 n (x,y) = m- 2n [d: Xo } z =o (C13) 
so that each of these terms corresponds to a superpotential deformation of the form [52] 



- I 9 n STr {$™ y F A F) (C.14) 
J s 

where the 4-cycle S is locally defined by the equation {z = 0}. Notice that Xo ^ 
log (h/ho)\s, and so the case 9q ^ corresponds to 

P yz \s = edJ p* x \ s = ed y 9 P*«\ s = (C.15) 

Now, in order to obtain the D-term one needs to compute 

P[e 2A Im *i] s A e F (C.16) 
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where P[a]s is the non-Abelian pull-back of a p-from a into the 4-cycle S. This means in 
particular that if a has two legs transverse to S, we should replace them with [& xy , $ xy ] 
(see section 2 of [53] for more details). Using this prescription, it is easy to see that the 
first two terms of the rhs of (C.10) give the usual D-term (2.2), with u = P[J]s- The third 
term of (C.10) would in principle give a correction to such D-term expression, however 
one can check that it is of order C(e 2 ). 16 

To summarize, from the results of [42] one can conclude that the expression for the 
D-term (2.2) remains unchanged after the non-perturbative effect has been taken into 
account. As shown in the main text this does not mean that the fluctuations satisfy 
exactly the same D-term equation before and after the non-perturbative effect. Indeed, 
while the differential equation is given by (4.23) in both cases, the F-term equations shift 
the values of ($), (A), and their conjugates, and so the actual D-term equations to be 
satisfied by the wavefunctions change. 

D Non-perturbative Yukawas and residues 

As pointed out in [18], and proved more generally in [54], one may compute the holomor- 
phic piece of the Yukawa couplings, that is the one that appears in the 4d superpotential, 
by means of an elegant residue formula. The purpose of this section is to explain how 
the analysis of [18] can be generalized to derive a residue formula for the general super- 
potential (4.1). We will focus on the case where only the holomorphic functions 6 and 
2 do not vanish, and show that the resulting residue formula gives Yukawa couplings in 
agreement with those obtained by directly computing wavefunctions overlaps, as done in 
the main text. 

D.l Non-perturbative F-terms and Yukawas 

One important result for deriving the residue formula for Yukawa couplings is to realize 

that the F-term equations are invariant under complexified gauge transformations, and 

that the holomorphic piece of the Yukawa couplings should be invariant under them [18]. 

16 Indeed, if we apply the prescription (C.16) to this (1,3) + (3, l)-form we always get an expression 
that contains either F ( 2, °) of D XjV & , which vanish at th order in e. 
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In the following we will generalize these concepts for the more involved superpotential 
(4.1), setting the stage for the residue formula discussed below. For simplicity we will 
focus in the case where only one term 6 n in (4.1) is non- vanishing, discussing in detail the 
cases n = and n — 1, but similar results can be derived for other values of n. 17 

The case 9 

Let us consider the superpotential 

W = W + W 1 = [ Tr($AF) + - I 9 Ti(F A F) (D.l) 
Js 2 J s 

that is, eq.(4.1) with only 9 ^ 0. From here we get the following set of equations of 
motion 

F(°' 2 ) = (D.2) 
d A <!> + ed6 A F {1 ' 1] = (D.3) 

that should be complemented with the D-term equation (4.23). This yields the following 
BPS equations for the background 

(F (0 ' 2) ) = (D.4) 
d {A) ($) + ed9 A(F^) = (D.5) 

and, expanding to first order in the fluctuation fields (</?, a), the zero mode equations 

d {A) a = (D.6) 
^^-^K^l + e^A^a + ^at) = (D.7) 

while the D-term gives 

u A (d {A) a + 8 {A) J) - l - ([($), <p] + [<p\ ($)]) = (D.8) 

17 As a byproduct of our computations, we will show that the the only fluctuations that enter into the 
Yukawa couplings are (a, ip), as opposed to their complex conjugates (a' , (p'). This assumption, taken all 
over the main text and in [22], is obviously true for the tree level superpotential (2.1), but not for the 
non-perturbative one (4.1), as the fluctuations a' appear explicitly in the latter. 
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One can check that the equations (D.6), (D.7) and (D.8) are invariant under the infinites- 
imal gauge transformations 

a -> a + B {A)X (D.9) 
a f -)• at + d {A)X (D.10) 
y? ¥>-i[<*),x] (D.ll) 
^ ^ + i [^) >x ] (D .i2) 

where x = is a re& l g au g e transformation parameter. However, if one is only interested 
in the F-term equations (D.6) and (D.7), then x can be arbitrary. This means that the 
equations are invariant under the complexified gauge group, as can be directly seen from 
the superpotential (D.l). 

Following [18], one finds that the general solution of (D.6) is 

a = 8 {A) £ (D.13) 

where £ is a scalar in the adjoint representation of complexified algebra. Using this result 
as well as the equation (D.5) one can rewrite (D.7) as follows 

d {A) (<p + i[(<l>),{;}-ede A(J-d {A) 0) = (D.14) 

from which we get the general solution 

<p = h- + ed9 A (at - d {A) £) (D.15) 

where h is an adjoint (2, 0)-form such that d^h = 0. 

As explained in the main text, the Yukawa couplings arise when we expand (D.l) to 
cubic order in fluctuations, plug in the solution to the equations of motion and perform 
the integration. We have two contributions to the Yukawa couplings, Y = Y + Y 1 , coming 
from Wq and Wi, respectively. The piece coming from W is 



Y = -i / Tr(v?AaAa) (D.16) 
Js 

— —i Tr((h-i[(^)^]-ed9 Ad {A} ^)AaAa) - ie Tr (d9 A a 1 " A a A a) 
■Is Js 



>s Js 
where we have used (D.15). The contribution coming from W\ is 



2 



Yi = J / fl Tr(FAF) (3 ) (D.17) 



s 
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where the subscript (3) means we only take the cubic terms in fluctuations. One can see 
that this expression amounts to 

Yl = -ie / 9 d {A) (Tr (a 1 " A a A a)) - ie 9 d {A) (Tr (a A a f A a f )) (D.18) 
Js Js 

Since when acting on a gauge invariant object, the operators d(A) an d <9(a) reduce to d 

and d, respectively, we can drop the subscript (A) in the last expression. Using that 9 is 

holomorphic we then obtain 

Yi = -ie [ o d(Tr(a ] AaAa)) - ie [ 8 (9 Ti (a A a f A a f )) (D.19) 
Js Js 

The second term in (D.19) is a boundary term so it vanishes by compactness of S. Adding 

the first term to the last term of (D.17) we have 

-ie / Tr (d9 Aa^ AaAa) - ie / 9 d (Tr (a' 1 " A a A a)) = -ie / d (Tr (0 O a f A a A a)) 
Js J s Js 

(D.20) 

which also vanishes upon integration. As a result, a* does not enter into the expression 
for the Yukawa coupling, which can be expressed as 

Y = -i [ Tr ([h -i[{$),Z]-ed0 o A d {A) £] A d {A) £ A d {A )Z) (D.21) 
Js 

To sum up, the expression for the Yukawa couplings for the case 9 ^ in this case is 
given by 

Y = -i / Tr(v?AaAa) (D.22) 
Js 

where a is giving by (D.13) and if by 

<p = h-i[{$),S]-ed6 Ad {A) S (D.23) 

which amounts to (D.15) with the fluctuation a) set to zero. In fact, for the purposes of 
computing Yukawa couplings to order (9(e) one may forget about the presence of a) in 
the superpotential, which indeed has been the working assumption of the main text. 

The case 9\ 

As a slightly more involved case let us consider the superpotential 

W = W + W 1 = [ Tr ($ xy F) AdxAdy+ e - [ 9,Tr {<$> xy F A F) (D.24) 
Js 1 Js 
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that corresponds to the case analyzed in 
equations of motion 18 



As in there, we have the following set of 



FAdxAdy + eOFAF = 
B A ^ + e[{F^\d A {9^ xy )} + {F^ 1 \d A {9^ xy )}] = 
{F^ 2 \d A {e^ xy )} + {F^\d A {e^ xy )} = 



(D.25) 
(D.26) 
(D.27) 



We now expand both 7-brane fields $ and A in the parameter e, as in (4.7), in order to 
write down the equations of motion for the fluctuations ip$, and a)^ 1 as 



(o) 



d {A) a {1) AdxAdy- i[(A) {1 \ a (0) ] A dx A dy + 6 { (F), d {A) a ] + d {A) a) 
d {A) ^-i[a^\ (<t>)(°)] -i[(A)^,^} + {{d {A) a + d {A) a\d {A) (9(<5> xy ))} (D.28) 



+{{F),d {A) {0<$> xy )-i[a\0{F)\}) 



(0) 







where we have made use of the zeroth order equations for the background. Following 
one can rewrite the above equations as 

* -X 1 ) .rwn^-HO) 



&aW = -6(K + iA) (o) 0^ - S7^ (0) 



(D.29) 



where ip is defined as in (5.1), and 6, K and A are defined by eqs.(3.19), (3.22) and 
(4.10) of [22]. Finally we have that S and p correspond to 

f {F°- z ,D y .}-6{F yy ,D r } 6{F xy ,D- z .}-{F°- z ,D y -} 0^ 
S = 9{F yx , D r } - {F yS , D x -} {F°- z , D x -} - 9{F XX , D- z -} 
y 0{F yy , D x -} — 9{F yx , D y -} 9{F xS , Dy-} — 9{F X y, D x -} J 



„t 



V J 



(D.30) 

Notice that eq.(D.29) corresponds to the zero mode equation (4.8), (4.9) of [22], except 
for the extra term proportional to ~~p* . As we will see in the following, in the computation 
of Yukawa couplings all terms including a) can be ignored, just like for the 9 case. Hence, 
in the equation of motion above we can ignore p, and so the zero mode analysis reduces 
to the one in [22]. 



Indeed, let us assume that we had found a wavefunction vector A such that 



->(o) ->(o) 
A = ip 



- -Hi) ~ ~ -Ho) 

d A A = -9(K + iA)6V 



(D.31) 



18 



In our conventions {A, B} = j{AB + BA). 
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which corresponds to solving the zero mode equation in 
with 



Since we can write S = (9a H 



H 



' 9{F yx ,-} —9{F xS , ■} ^ 

^ { ^yy > '} — 9{F xy , ■} 

v {f;- z ,-} -{f°»-} o j 

we have that a general solution to (D.29) is given by 



(D.32) 



(D.33) 



which is analogous to the solution (D.15) for the case Oo, and we can basically repeat the 
steps of this case to show that ~~ff only appears in the trilinear couplings through total 
derivatives. On the one hand, we have a trilinear term involving ~ff that arises from Wo 

Y Die j Tr(y? A [w, a] + {a\d {A )(9(^ xy ))} A a A a) (0) (D.34) 
Js 

where we have defined vo = (Hp)idx + (Hp) 2 dy = {a\ F^ 1 ' 1 ^}^. On the other hand, W\ 
contains the following terms involving a* 



-it 



Tr [9 <p xy ({(F), [at, a] } + % [d {A) a, 5 {A) <f] + id {A) <J A d {A) J) (D.35) 



+ 6{$ xy ) {{d {A) a\ a A a} + {d {A) a, [a\ a] } + {d {A) a\ [a\ a] }) ) 
Putting all these terms together one finally obtains 



(0) 



Tr (d ([a , 0($ xy )a i ]Aa i + id^ xy (a f A do) + [a\ da})) + 2d {{a\ 9(<f> xy ) } A a A a)) 



(0) 



(D.36) 

which vanishes upon integration. Therefore, one can set a) = from the very beginning 
without loss of generality, and recover the analysis performed in [22]. 



D.2 Residue formula 

In this section we will derive the residue formula for the holomorphic Yukawas induced by 
the non-perturbative terms in the general superpotential (4.1). We will treat separately 
the dependence on 9 and 9 2 , which are the parameters of interest for the S'0(12) model 
described in section 3. 
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Yukawas due to 9 

In section D.l we have seen that the Yukawa couplings including corrections are given by 

Y = —im* / Tr (ip A a A a) 
Js 

= m*fabc / [fc A a a A a b + ip b A a c A a a + (p a A a b A a c ] (D.37) 
Js 

These couplings can be conveniently rewritten by using the F-terms found after shifting 
the background as explained in section 4.2. In particular, eq. (4.20), see also eq. (5.5), 
tells us that 

dif p = -im 2 q^{p)a p dx Ady- ed9 A da p + 0(e 2 ) (D.38) 
On the other hand, eq. (4.19) amounts to 

a p = dt; p + 0(e 2 ) (D.39) 

where £ p is a function with properties to be discussed shortly. Clearly (D.38) can then be 
integrated to obtain 

(f p = hpdx A dy — im 2 qq>(p)^pdx A dy — ed9o A d£ p (D.40) 

where h p is a holomorphic function. 

Substituting the F-terms in the coupling Y we find 

Y = m*f abc {h c a a A a b A dx A dy + 8 [((p b A a a - ip a A a b )£ c ] - ed [9 d(a a A a b £ c )]} 
Js 

(D.41) 

to 0(e 2 ). Upon integration the total derivatives lead to boundary terms that vanish 
because, as shown in section 5.1, the zero modes for a p and ip p are localized on the curve 
E p . Taking into account that h c is holomorphic, and using (D.39) once more, then gives 

Y = mj abc [ 8 (hctadtb) A dx A dy (D.42) 
Js 

Before further integrating we will first study the functions £ p in more detail. 

The idea is to determine £ p recursively from (D.40). We know that ip p has an expansion 
(5.7) in powers of e, and in fact the corrections ip^ are explicitly given in section 5.1. It 
is natural to expand £ p in a similar way, namely 

Zp = + e£« + 0(e 2 ) (D.43) 
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From (D.40) it then follows that 

^ =-^T-M ] -K) (D.44a) 

= W + dJoUf - d y 9 Q d x if] (D.44b) 

The functions £ p are requested to be regular at the locus q$(p) = 0. Recall that q<s>(b + ) = y 
whereas q<p(a + ) = —x. At zeroth order, to avoid poles it must be that h p = ^p - 1 1^=0- For 
example, h a = m*fi{y) and h b = m*gj{x). We will later study the regularity of the 
corrections 

To arrive at the residue formula it is important to realize that after integrating the 
total derivative in (D.42) the terms involving ip p in £ p will not contribute to the couplings. 
The justification to drop these terms is again the localization of ip p on the curve S p . It 
is then useful to define auxiliary functions r\ p such that £ p reduces to — r\ p away from the 
curve, as it is done in [18] where i] a and r] b are denoted £12 and (23 respectively In practice 
i] p is found discarding the ip p dependence in (D.44). Therefore, 



% . e 
= 2 — T^ h p + 



m 2 q^(p) m 4 g$(p) 



d x o d y I — j— - J — d y 9 d x ^ 



C(e 2 ) (D.45) 



Notice that rj p is holomorphic. Thus, after integrating the total derivative in (D.42) the 
coupling reduces to 

Y = mj abc / 8 (h c r] a ^ b ) AdxAdy (D.46) 

JdS 

Integrating the last total derivative then gives 

Y = m*f abc / (h c r] a r] b ) dx A dy (D.47) 
J-r 

where the final integration region is effectively the product of two circles, \x\ = Si and 
\y\ = S2, each of infinitesimal radius. Evaluating by residues gives 19 

Y = mJ abc 7i 2 Res(h c r] a r] b ,x = 0,y = 0) (D.48) 

This is the residue formula. The couplings only depend on the holomorphic functions 
h p that determine r\ p through eq. (D.45). We will shortly present explicit examples of 
couplings computed via residues. 



19 We also divide by (2i) 2 to match the normalization of dvols used in evaluating overlap integrals of 
zero modes. 
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Regularity of ^ 

The ^p 1 "* are easily found inserting the corrected wavef unctions in (D.44). We will 
briefly discuss the sectors a + , which is analogous to b + , and c + . To keep track of the 
generations in the a + and b + sectors we add superscripts i and j. To connect with section 
5.1 we take 9 = im 2 (9 Q + x9 0x + y0 0y ). 

Lets us consider the a + sector. The correction to the wavefunction (p l a+ is given in 
(5.9). Substituting in (D.44), together with (p % Jp = m*e XaXUa fi(v a ) and h l a+ = m*/j(y), 
and Taylor expanding, leads to 



m 2 x 



A i (y) + ^(29 0x -C a e 0y )f:\y) 



+ ■■■ (D.49) 



where the ellipsis stands for terms regular at x — 0. Hence, to avoid the pole at x — it 
suffices to take 

Mv) = | (Caeoy- 29 0x ) fl'(y) (D.50) 

up to a holomorphic function that vanishes at x — 0. Since fi(y) = m£~^ ^\y^ z ~ l \ we see 
that only Ai ^ and moreover 

M = 7^a ; a = mlC a (C, a 9 Qy - 29 0x ) (D.51) 

as reported in section 5.1. The results for the b + sector are analogous. For to be free 
of poles B 2 = B 3 = 0, whereas B\ = 7^60, with b = ml(b((b9o x — 29 0y )- 

We now address the c + sector. The corrected wavefunction is given in (5.20). In this 
case iff} = m^ c e ( ' c( - x ~ y ^ Uc , so that h c + = m^ c . Inserting in (D.44) we arrive at 

& = ~2T S m *^ C + ■■■ ( D - 52 ) 

m z [x — y) 

where • • • represents terms regular at x — y. Hence, to remove the pole at x = y it is 
enough to set C = 0, as we announced in section 5.1. 

Couplings 

We have just seen that h c + = m^ c . Then, according to (D.48), to compute the 
couplings Y lj to 0(e 2 ) we just need to extract the residue of rj l a + r li+ & t {x,y) = (0,0). 
Recall that the a + and b + sectors are analogous. We will present two illustrative examples. 
Other couplings can be found in a similar way. 
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Let us first calculate Y 33 . For the third generation we have h 3 a+ = m*7 3 , and likewise 
h 3 + = m^f. From (D.45) we find 

* = 4- M + f!h ^) + < D53 > 

rrrrr \ x 2 J 

* = -^(m^ + ^I^W*) (D.54) 
Extracting Res(r/^+ 7/ 3 +, x = 0, y = 0) gives 

Y$ b+C+ = tt 2 (^-Yf a+b+c+ lcl hl + 0( e 2 ) (D.55) 

This coupling remains then uncorrected to 0(e 2 ), as we also found calculating the overlap 
integral of corrected wavef unctions. 

We now compute Y 13 . Recall that f\ = ml'jKy + ( a x) 2 , so that h l a+ {y) = ml^y 2 . 
Inserting in (D.45) yields 



X \ X 



m^x 

1 ^3 



+ C(e 2 ) (D.56) 



r] 3 + is still given by (D.54). Evaluating the residue of i] a+ ri 3 + at (0,0) we find that 



f I T) 

y^ +c+ = -^* 2 fa +b+c+ IdhiOo* + 0(e 2 ) (D.57) 



6 



in agreement with the coupling (6.11c) obtained from the overlap integral of corrected 
zero modes. 

Yukawas due to 9 2 

We now start from the superpotential 

W = W + W 1 = [tt (<S> xy F) AdxAdy + ^ [ 6 2 STr ($ 2 F A F) (D.58) 
Js 2 Js 

We will take 6 2 constant. The non-perturbative piece W\ includes a cubic term, shown 

in (6.4), depending on the fluctuations ((p,a). In turn this term implies the corrected 

Yukawa coupling (6.9) which can be written as Y\ = Yi<$) + Yi{f)- In the 50(12) model, 

taking into account the group properties (6.14), Y ± ^ is given in simplified notation by 

Y 1( $) = % -em*$ abc / 6 2 [Da a A Da b Lp cxy + Da c A Da a (p bxy + Da b A Da c ip axy ] (D.59) 
1 Js 
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where 6 2 is the coordinate dependent parameter 

§ 2 = ^m 2 (y-2x) (D.60) 
Similarly, for (F) = 3Qf, the flux piece Y^ F ) reduces to 

Yl(F) = ^mj abc J 9 2 ((PaxyPbxyDac + <P cxy<P axyD d h + (PbxyVcxyDaa) A # (D.61) 

The total coupling can be expanded in powers of e as 

Y = Y + Y m + Y l[F) + 0(e 2 ) (D.62) 

where Y arises from the tree-level superpotential Wo and it is actually given in (D.37). 

In the n = case discussed before we have seen that to arrive at a residue formula we 
can use the F-terms to expose the total derivatives in the integrand in the full coupling 
Y. However, in principle this approach is now more difficult because the F-terms (4.21) 
and (4.22) are formidable equations to work with. A way to overcome the difficulties is 
to switch to variables A m and & xy by performing the generalized SW map introduced in 
(4.11). We will then first take a detour to show that in the new variables the F-term 
equations, as well as the Yukawa couplings, are much simpler and in fact independent of 
fluxes. 

With 6> 2 constant, and other 9 n = 0, the generalized SW map (4.11) is given by 

Afh = An + = A fh + e6 2 S {<Z> xy [A x (d y A m + Fyjn) - AyidxArn + F xfh )}} (D.63a) 

K = ®xy + $xy = $xy + fa [A x (d y + D y ){^ xy ) - A y (d x + D X )(&J\ (D.63b) 

up to 0(e 2 ) corrections. We also have to expand in fluctuations ((p xy , %) defined by 

®xy = (&xy) + 0xy 5 = (Afh) + % (D.64) 

In turn the vevs and the fluctuations have expansions in powers of e of the form 

(K) = (**y) {0) + (K) ll) + O(<?) (D.65) 
In the holomorphic gauge (Afn)^ = the vevs reduce to 

($xy) = (<5>xy) {0) + O(e 2 ) ; (A rh ) = + O(e 2 ) (D.66) 
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as can be proved using (4.12) and (4.13). For the fluctuations we find 

0xy = ¥xy + <fxy ] «m = «m + (D.67) 

where 

xy/^Pxy ) - (A y )d x (($ xy )<p xy )} (D.68a) 
a„ = 2e6 2 S [($ xy ) {(A x )d y a m - (A y )d x a m )} 

+ e6 2 S [p xy {(A x )(F m ) - (A y ) (F xfh ))] (D.68b) 

up to 0(e 2 ) terms. 

To derive the equations of motion satisfied by (p xy and d„ we substitute (D.67), to- 
gether with (D.68), into (4.21) and (4.22). In the end we obtain 

dx(i y — d y d x = (D.69a) 
d^ xy + t[(^ xy ) (0) ,a M } = 2e9 2 S [{® xy ) (d x a r?l d y (^ xy ) - d^d^ xy ))f ] (D.69b) 

up to 0(e 2 ) terms. Notice that these equations are independent of fluxes. In fact they can 
be derived from (4.21) and (4.22) ignoring flux terms and replacing (a™, f xy ) by (dm, <fixy)- 
For instance, using (3.8) from (D.69b) we deduce 

d m <p pxy = -im 2 q^(p)d pm - 2e9 2 m 4 [t x (p)d y d pm - t y (p)d x d pm ] + 0(e 2 ) (D.70) 

where the t m (p) are defined in (5.39). This equation corresponds to the F-term (5.38), 
which together with the equivalent of (D.69a), was applied in section 5.3 to find the zero 
modes in the absence of fluxes. Thus, the wavefunctions for ip pxy and d pm can be inferred 
from the results in section 5.3. In particular, we see that these wavefunctions are localized 
on the curves S p . 

We now continue as in the n = case. To begin notice that (D.70) can be written as 

d(p p = -im 2 q^(p)a p - 2e0 2 m A t(p) A da p + £>(e 2 ) (D.71) 

Recall also that dt(p) = 0. Since (D.69a) means a p = d^ p + 0(e 2 ) we can integrate (D.71) 
to obtain 

ip p = h p dx A dy - im 2 q^(p)i p dx Ady- 2e9 2 m 4 t(p) A di p + 0{t 2 ) (D.72) 
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where h p is a holomorphic function. We can solve for £ p iteratively and then define 
auxiliary functions fj p such that £ p reduces to —fj p on the curve S p . It follows that 



-h 



2e9 9 



m 2 q$(p) q$(p) 



tx(p)dy 



h, 



- ty{p)d x 



h, 



0(e 



(D.73) 



" \QM J " yyr '~ x \ q$(p) 
We see that the fj p are also holomorphic. 

Let us now discuss the Yukawa couplings. To begin we will express them in terms of 
the hat variables by substituting (D.67) in the full expansion (D.62). The result hinges 
crucially on the corrections (D.68) in which we can insert the vevs to order e°. In the 
50(12) model, with (A m ) = 2l m Q F , we find 



<fp Xy = -2em 2 9 2 [% x d y (r p (p pxy ) - %d x (r p ip 



pxyj 



(D.74) 



Here r p = xs xx (p) + (x + y)s xy (p) + ys yy (p), where the s mn (p) are provided in table 1. 
Similarly 

a P rn = -2e9 2 m 2 r p [% x d y a pr - n - 21^0^] - e0 2 s p [% x drn% - %drn% x ] ip pxy (D.75) 

where s p = s xx (p) + 2s xy (p) + s yy (p). Since we are working in the holomorphic gauge we 
have written ^ nfh = —dm^i n for the flux components. 

After long but direct calculations we find that the full Yukawa couplings take the form 



Y = Y + Y m + 0(e' 



(D.76) 



where Y and Y^ are obtained by replacing (a p , ip p ) with (a p , <p p ) in the formulas (D.37) 
for Y , and (D.59) for Yi^y Moreover, in Y^ all flux dependence disappears because 
only normal derivatives da p do enter in the final expression. In deriving (D.76) we have 
used the F-terms (D.69). We have also dropped boundary terms that vanish due to 
localization of the wavef unctions. 

The next step is to use repeatedly the equations satisfied by a p and (p p to show that 
(D.76) implies 

ie§ 2 



Y = TO*/, 



abc I 

Js 



a a A db A dx A dy H dd a A ddb 



(D.77) 



In the proof we have again discarded boundary terms. The result also relies on the 
properties 

d{q„(b + )6 2 ) = Am 2 9 2 [t{a + ) - t(c + )] ; d{q^{a + )9 2 ) = 4m 2 6 2 [t{c + ) - t(b + )] (D.78) 
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which are valid in the 50(12) model. At this point we can use d p = d£ p + 0{e 2 ) to 
continue integrating. Furthermore, at the boundary we can substitute £ p by fj p thanks to 
localization. In this way we find 



Y = m*/ ( 



abc 



h, 



n 



VaVb 



dx A dy 



(D.79) 



Evaluating by residues as in the n = case then gives 



Y = m*f abc ir' z Res (k, 



VaVb - l -^r {dxfjadyfn, - d y fj a d x 7] b ) 



,x = 0,y = 



(D.80) 



This is our final residue formula. 

Couplings can now be easily calculated. To this purpose we first determine fj p from 
(D.73). The holomorphic functions h p must be such that ^ is regular at q$(p) = 0. This 
implies 



;(o) 



g$(p)=0 



(D.81) 



hp ~ ¥ 'pxy\q<s>(p)=0 — { P pX y\ 

Therefore, h a +{y) = m*fi(y), h b +(x) = m*gj(x), and h c + = m*7 c . From the results in 
section 5.3 we can further show that the resulting ffi are indeed regular. 

Since h c + = m^ c , to compute the couplings Y lj to 0(e 2 ) it suffices to determine 
fj l a+ and ff b+ . As a first example let us consider Y 33 . Substituting h 3 a+ (y) = m*7^, and 
h 3 + (x) = m^l in (D.73) yields 



m?x 



2x 3 



■(x-y) 



m 2 y 



(D.82) 



to 0(e 2 ). Plugging in (D.80) readily reproduces the uncorrected coupling in (D.55). For 
a more interesting example we calculate Y 32 . With h 2 + = ml^x we find 

-ix ed 2 



"2 2 2 

Vb+ = m *% 



m 2 y 6y 



(D.83) 



From (D.80) we now obtain 



V-32 

*a+b+c+ 



iem^ii 2 



Qm 2 



f a+ b + c + icilile 2 + o{t 2 ) 



(D.84) 



which matches (6.15b) taking 9 2 = i9 20 . It is easy to check that the coupling Y 23 derived 
from the residue formula agrees with (6.15a), and that all other couplings vanish to 0(e 2 ). 



92 



References 



[1] For an overview of string phenomenology see L.E. Ibanez and A.M. Uranga, String 
Theory and Particle Physics. An Introduction to String Phenomenology, Cambridge 
University Press (2012). 

[2] G. Aldazabal, L. E. Ibanez, F. Quevedo and A. M. Uranga, "D-branes at singularities: 
A Bottom up approach to the string embedding of the standard model, " JHEP 0008, 
002 (2000) [arXiv:hep-th/0005067]. 

[3] R. Donagi, M. Wijnholt, "Model Building with F-Theory, " [arXiv:0802.2969 [hep-th]]. 

[4] C. Beasley, J. J. Heckman and C. Vafa, "GUTs and Exceptional Branes in F-theory 

- I," JHEP 0901 (2009) 058 [arXiv:0802.3391 [hep-th]]. 

[5] C. Beasley, J. J. Heckman and C. Vafa, "GUTs and Exceptional Branes in F-theory 

- II: Experimental Predictions," JHEP 0901 (2009) 059 [arXiv:0806.0102 [hep-th]]. 

[6] R. Donagi and M. Wijnholt, "Breaking GUT Groups in F-Theory, " [ arXiv:0808.2223 
[hep-th]]. 

[7] J. J. Heckman, "Particle Physics Implications of F-theory," arXiv: 1001.0577 [hep- 
th]; T. Weigand, "lectures on F-theory compactifications and model building," 
Class. Quant. Grav. 27, 214004 (2010) [arXiv: 1009.3497 [hep-th]]; L. E. Ibanez, 
"From Strings to the IHC: les Houches lectures on String Phenomenology, " 
arXiv:1204.5296 [hep-th]; G. K. Leontaris, "Aspects of F-Theory GUTs," PoS 
CORFU 2011 (2011) 095 [arXiv: 1203.6277 [hep-th]]. 

[8] H. Hayashi, T. Kawano, R. Tatar and T. Watari, "Codimension-3 Singularities and 
Yukawa Couplings in F-theory," Nucl. Phys. B 823, 47 (2009) [arXiv: 090 1.4941 [hep- 
th]]. 

[9] H. Hayashi, T. Kawano, Y. Tsuchiya and T. Watari, "Flavor Structure in F-theory 
Compactifications," JHEP 1008, 036 (2010) [arXiv:0910.2762 [hep-th]]. 

[10] J. J. Heckman and C. Vafa, "Flavor Hierarchy From F-theory," Nucl. Phys. B 837 
(2010) 137 [arXiv:0811.2417 [hep-th]]. 

93 



[11] A. Font and L. E. Ibanez, "Yukawa Structure from U(l) Fluxes in F-theory Grand 
Unification," JHEP 0902, 016 (2009) [arXiv:0811.2157 [hep-th]]. 

[12] E. Dudas and E. Palti, "Frog gatt- Niels en models from E(8) in F-theory GUTs," 
JHEP 1001, 127 (2010) [arXiv:0912.0853 [hep-th]]. 

[13] G. K. Leontaris and G. G. Ross, "Yukawa couplings and fermion mass structure in 
F-theory GUTs, " JHEP 1102, 108 (2011) [arXiv: 1009.6000 [hep-th]]. 

[14] S. Krippendorf, M. J. Dolan, A. Maharana and F. Quevedo, "D-branes at Toric 
Singularities: Model Building, Yukawa Couplings and Flavour Physics, " JHEP 1006, 
092 (2010) [arXiv: 1002. 1790 [hep-th]]. 

[15] D. Cremades, L. E. Ibanez and F. Marchesano, "Yukawa couplings in intersecting 
D-brane models," JHEP 0307, 038 (2003) [arXiv:hep-th/0302105]. 

[16] D. Cremades, L. E. Ibanez and F. Marchesano, "Computing Yukawa couplings from 
magnetized extra dimensions," JHEP 04 05, 079 (2004) [arXiv:hep-th/0404229]. 

[17] F. Marchesano and G. Shiu, "MS SM vacua from flux compactifications," Phys. Rev. 
D 71, 011701 (2005) [arXiv:hep-th/0408059]. "Building MS SM flux vacua," JHEP 
0411, 041 (2004) [arXiv:hep-th/0409132]. 

[18] S. Cecotti, M. C. N. Cheng, J. J. Heckman and C. Vafa, "Yukawa Couplings in 
F-theory and Non-Commutative Geometry," [arXiv:0910.0477 [hep-th]]. 

[19] J. P. Conlon and E. Palti, "Aspects of Flavour and Supersymmetry in F-theory 
GUTs," JHEP 1001, 029 (2010) [arXiv:0910.2413 [hep-th]]. 

[20] A. Font and L. E. Ibanez, "Matter wave functions and Yukawa couplings in F-theory 
Grand Unification," JEEP 0909, 036 (2009) [arXiv:0907.4895 [hep-th]]. 

[21] F. Marchesano and L. Martucci, "Non-perturbative effects on seven-brane Yukawa 
couplings," Phys. Rev. Lett. 104, 231601 (2010) [arXiv:0910.5496 [hep-th]]. 

[22] L. Aparicio, A. Font, L. E. Ibanez, F. Marchesano, "Flux and Instanton Effects in 
Local F-theory Models and Hierarchical Fermion Masses," JHEP 1108 (2011) 152 
[arXiv: 1104.2609 [hep-th]]. 



[23] A. Collinucci, "New F-theory lifts. II. Permutation orientifolds and enhanced singu- 
larities," JHEP 1004, 076 (2010) [arXiv:0906.0003 [hep-th]]. 

[24] R. Blumenhagen, T. W. Grimm, B. Jurke and T. Weigand, "F-theory uplifts and 
GUTs," JHEP 0909, 053 (2009) [arXiv:0906.0013 [hep-th]]. "Global F-theory GUTs," 
Nucl. Phys. B 829, 325 (2010) [arXiv:0908.1784 [hep-th]]. 

[25] J. Marsano, N. Saulina and S. Schafer-Nameki, "Monodromies, Fluxes, and Compact 
Three- Generation F-theory GUTs," JHEP 0908, 046 (2009) [arXiv:0906.4672 [hep- 
th]]. 

[26] C. Cordova, "Decoupling Gravity in F-Theory," Adv. Theor. Math. Phys. 15, 689 
(2011) [arXiv:0910.2955 [hep-th]]. 

[27] J. Marsano, N. Saulina and S. Schafer-Nameki, "Compact F-theory GUTs with U(l) 
(PQ)," JHEP 1004, 095 (2010) [arXiv:0912.0272 [hep-th]]. 

[28] T. W. Grimm, S. Krause and T. Weigand, "F-Theory GUT Vacua on Compact 
Calabi-Yau Fourfolds," JHEP 1007, 037 (2010) [arXiv:0912.3524 [hep-th]]. 

[29] T. W. Grimm, "The N=l effective action of F-theory compactifications, " Nucl. Phys. 
B 845, 48 (2011) [arXiv: 1008.4133 [hep-th]]. 

[30] E. Dudas and E. Palti, "On hypercharge flux and exotics in F-theory GUTs," JHEP 
1009, 013 (2010) [arXiv: 1007. 1297 [hep-ph]]. 

[31] H. Abe, T. Kobayashi and H. Ohki, "Magnetized orbifold models," JHEP 0809, 043 
(2008) [arXiv: 0806. 4748 [hep-th]]; 

H. Abe, K. S. Choi, T. Kobayashi and H. Ohki, "Three generation magnetized orbifold 
models," Nucl. Phys. B 814, 265 (2009) [arXiv:0812.3534 [hep-th]]. 

[32] J. P. Conlon, A. Maharana and F. Quevedo, "Wave Functions and Yukawa Couplings 
in Local String Compactifications," JHEP 0809, 104 (2008) [arXiv:0807.0789 [hep- 
th]]. 



95 



[33] P. Di Vecchia, A. Liccardo, R. Marotta, F. Pezzella, "Kahler Metrics and Yukawa 
Couplings in Magnetized Brane Models," JHEP 0903 (2009) 029. [arXiv:0810.5509 
[hep-th]]. 

[34] P. G. Camara and F. Marchesano, "Open string wavef unctions in flux compactifica- 
tions," JHEP 0910, 017 (2009) [arXiv:0906.3033 [hep-th]]. "Physics from open string 
wavefunctions," PoS E PS-HEP2009, 390 (2009). 

[35] H. Georgi, Lie algebras in Particle Physics: from isospin to unified theories, Westview 
Press, 1999. 

[36] E. Palti, "Wavefunctions and the Point of E$ in F-theory," JHEP 1207, 065 (2012) 
[arXiv: 1203.4490 [hep-th]]. 

[37] L. E. Ibahez, F. Marchesano, D. Regalado and I. Valenzuela, "The Intermediate 
Scale MSSM, the Higgs Mass and F-theory Unification," JHEP 1207, 195 (2012) 
[arXiv: 1206.2655 [hep-ph]]. 

[38] D. Baumann, A. Dymarsky, I. R. Klebanov, J. M. Maldacena, L. P. McAllister and 
A. Murugan, "On D3-brane Potentials in Compactifications with Fluxes and Wrapped 
D-branes," JHEP 0611, 031 (2006) [hep-th/0607050]. 

[39] S. A. Abel and M. D. Goodsell, "Realistic Yukawa couplings through instantons in 
intersecting brane worlds," JHEP 0710, 034 (2007) [arXiv:hep-th/0612110]. 

[40] M. Bianchi, A. Collinucci and L. Martucci, "Magnetized E3-brane instantons in F- 
theory," JHEP 1112, 045 (2011) [arXiv: 1107.3732 [hep-th]]. 

[41] M. Cvetic, I. Garcfa-Etxebarria and J. Halverson, "Three Looks at Instantons in F- 
theory - New Insights from Anomaly Inflow, String Junctions and Heterotic Duality, " 
JHEP 1111, 101 (2011) [arXiv: 1107.2388 [hep-th]]. 

[42] A. Dymarsky and L. Martucci, "D-brane non-perturbative effects and geometric de- 
formations," arXiv: 1012.4018 [hep-th]. 

[43] N. Seiberg and E. Witten, "String theory and noncommutative geometry, " JHEP 
9909, 032 (1999) [arXiv:hep-th/9908142]. 



[44] G. Ross and M. Serna, "Unification and fermion mass structure," Phys. Lett. B 664 
(2008) 97 [arXiv:0704.1248 [hep-ph]]. 

[45] G. Elor, L. J. Hall, D. Pinner and J. T. Ruderman, "Yukawa Unification and the 
Superpartner Mass Scale," JHEP 1210 (2012) 111 [arXiv: 1206.5301 [hep-ph]]. 

[46] H. Georgi, C. Jarlskog, "A New Lepton - Quark Mass Relation in a Unified Theory, " 
Phys. Lett. B86 (1979) 297-300. 

[47] K. Bora, "Updated values of running quark and lepton masses at GUT scale in SM, 
2HDM and MSSM," arXiv: 1206.5909 [hep-ph]. 

[48] P. G. Camara, E. Dudas and E. Palti, "Massive wavefunctions, proton decay and 
FCNCs in local F-theory GUTs," JHEP 1112, 112 (2011) [arXiv: 11 10.2206 [hep- 
th]]. 

[49] D. Baumann, A. Dymarsky, S. Kachru, I. R. Klebanov and L. McAllister, "D3-brane 
Potentials from Fluxes in AdS/CFT, " JHEP 1006, 072 (2010) [arXiv: 100 1.5028 [hep- 
th]]. 

[50] R. C. Myers, "Dielectric branes, " JHEP 9912, 022 (1999) [arXiv:hep-th/9910053]. 

[51] P. Koerber and L. Martucci, "From ten to four and back again: how to generalize 
the geometry," JHEP 0708, 059 (2007) [arXiv:0707.1038 [hep-th]]. 

[52] L. Martucci, "D-branes on general N=l backgrounds: Superpotentials and D -terms," 
JHEP 0606, 033 (2006) [hep-th/0602129]. 

[53] F. Marchesano, P. McGuirk and G. Shiu, "Chiral matter wavefunctions in warped 
compactifications," JHEP 1105, 090 (2011) [arXiv:1012.2759 [hep-th]]. 

[54] S. Cecotti, C. Cordova, J. J. Heckman and C. Vafa, "T-Branes and monodromy," 
JHEP 1107, 030 (2011) [arXiv: 1010.5780 [hep-th]]. 



97 



